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Abstract 



Structures of light unstable nuclei, Li, Be, B, and C isotopes are systematically 
studied with a microscopic method of antisymmetrized molecular dynamics. 
The theoretical method is found to be very useful to study ground and excited 
states of various nuclei covering unstable nuclei. The calculations succeed 
to reproduce many experimental data for nuclear structures; energies, radii, 
magnetic dipole moments, electric quadrupole moments, transition strength. 
In the theoretical results it is found that various exotic phenomena in unstable 
nuclei such as molecular-like structures, neutron skin, and large deformations 
may appear in unstabel nuclei. We investigate the structure change with the 
increase of neutron number and with the increase of the excitation energies, 
and find the drastic changes between shell-model-like structures and clustering 
structures. The mechanism of clustering developments in unstable nuclei are 
discussed. 
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I. INTRODUCTION 



Owing to the progress of the experimental technique, the information of 
the ground and excited states of unstable nuclei is increasing rapidly. There 
are various interesting subjects which are characteristics of unstable nuclei. 
Our aim is to systematically study the nuclear structures of light nuclei cov- 
ering stable and unstable regions in order to understand the features of the 
nuclear many-body system. In the light nuclear region, ground-state proper- 
ties have become known experimentally for the unstable nuclei up to the drip 
line. The experimental data give various kinds of information such as binding 
energies, radii, magnetic dipole moments, and electric quadrupole moments, 
and so on jj]-13|. By the help of these experimental data, many interest- 
ing phenomena of the structures of the unstable nuclei have been suggested; 
neutron halo and skin structures, vanishing of the magic number, abnormal 
spin-parity of the ground state, clustering structures, large deformations in 
unstable nuclei. 

In many of the theoretical studies for light nuclei, inert cores or clusters 
have been assumed. For example, three-body models have been applied to 
11 Li by regarding it as a 9 Li+2n system and also applied to 6 He as an a+2n 



system |14J-[n|. In these studies the neutron-halo structure has been theo- 
retically investigated. In the studies of very light nuclei with A < 10 with 
extended cluster models in which a or t cluster cores and surrounding nu- 
cleons are assumed [|l7| -|l9||, it is suggested that the clustering structure may 
appear also in unstable nuclei. However, they have not yet reached to the 
systematic investigations covering the wide region of nuclei, since it is diffi- 
cult to apply these models to the study of general heavier unstable nuclei. 
For heavier unstable nuclei, systematic studies have been done by using other 
theoretical frameworks such as shell model approaches [BO] and the methods 
of mean field theory [21,22]. Some of them suggested new features such as 



large deformations and neutron skin structures. However the applicability 
of mean-filed approaches is not necessarily assured in light nuclei because of 
possible clustering structures. 

It is already well known that clustering structures p3|-p7[| appear in the 
ground states of ordinary light nuclei with iV ~ Z as seen in the t + a clus- 
ter structure of 7 Li and in the a + a cluster structure of 8 Be, and in the 
16 0+q cluster structure of 20 Ne. Therefore, in light unstable nuclei one of 
the problems to be solved is the clustering structure. Is the clustering seen 
also in unstable nuclei ? If it is the case, what is the feature of the clusters. 
How does the developed clustering structure of stable nuclei change with the 
increase of the neutron number in a series of isotopes ? It is naturally ex- 
pected that the clustering structures of unstable nuclei are found not only 
in the ground state but also in the excited states. Recent experimental data 
suggest that the clustering structure may appear in excited states of unstable 
nuclei. For example, many excited states of 12 Be have been found with the 
experiments of the breakup reactions from 12 Be to 6 He+ 6 He and 8 He+ 4 He 
channels p8[|. In theoretical studies, many groups have suggested cluster- 
ing structures in neutron-rich nuclei and tried to discuss the features of the 
clustering [17-[lS,2S,3C]. It is essential for the systematic studies of unstable 



nuclei to describe both clustering aspects and shell-model-like aspects in one 
theoretical framework which is applicable to general nuclei. We are interested 
in how the structure changes with the increase of the neutron number and 
the excitation energy. 

Our aim is to make systematic investigations of the structures of the 
ground and excited states of unstable nuclei with a microscopic model which 
is free from such assumptions as the inert core or the existence of clusters. Re- 
marking the development of clustering structure, we try to understand many 
characteristic phenomena seen in unstable nuclei. We discuss the mechanism 
of clustering developments of unstable nuclei. 

In this paper, we adopt a theoretical method of antisymmetrized molec- 
ular dynamics(AMD). Ono et al. have developed the method of AMD for 



the study of nuclear reactions [31-35]. The framework of AMD has been 



extended by Kanada-En'yo (one of the author of this paper) et al., and has 
been applied to the studies of nuclear structures [p9]^0|j3^ -|39| . From these 
studies AMD has already proved to be a very useful theoretical approach for 



investigating the structure of the ground and excited states of light nuclei. In 
the AMD framework, basis wave functions of the system are given by Slater 
determinants where the spatial part of each single-particle wave function is 
a Gaussian wave packet. One of the characteristics of AMD is the flexibility 
of the wave function which can represent various clustering structures and 
shell-model-like structures, which is because no inert cores and no clusters 
are assumed. Another characteristic point of AMD is the frictional cooling 
method which is adopted in the energy variation for obtaining the ground 
and excited states. In the case of the simplest version of the AMD [29|,30| 



in this paper, the energy variation is made after the parity projection but 
before the total-spin projection. With this simplest method, we can describe 
low-lying levels of the lowest bands with positive and negative parities. For 
the study of the excited states we adopt the AMD approach which users the 
variation after the spin-parity projection (VAP). VAP calculations in AMD 
framework have been already found to be advantageous for the study of the 
excited states of light nuclei 3S.3S]. By the use of microscopic calculations 



with the obtained wave functions, we can easily acquire the theoretical values 
of energies, radii, magnetic dipole moments, electric quadrupole moments, 
and the transition strength of El, El, and (5 which are useful quantities to 
deduce the informations of structures from the experimental data. 

This paper is organized as follows. In the next section (Sec. ||) we explain 
the formulation of AMD for the study of nuclear structure. The effective 



interactions are described in Sec. |III| . In Sec. IV we show the results and 
give discussions of the Li, Be, B, and C isotopes based on the calculations 
of the simplest version of AMD. The study with VAP calculations in the 
framework of AMD is reported in Sec. [V|, where we discuss the structure of 



the excited states of neutron-rich Be isotopes. In Sec. VI, we mention about 



the mechanism of clustering developments of Be isotopes. Finally we give the 



summary in Sec. VII 



II. FORMULATION 

The AMD (antisymmetrized molecular dynamics) is a theory which is 
applicable to the studies of the nuclear structure and the nuclear reaction. 
Here we only explain the AMD framework for the study of nuclear structures. 
As for the AMD theory for the study of nuclear reaction, the reader is referred 



to Ref. [32]. 



A. AMD wave function 

In AMD framework, the wave function of a system is written by a linear 
combination of AMD wave functions, 

$ = c$amd + c&amd + ■■■■ (!) 



An AMD wavefunction is a Slater determinant of Gaussian wave packets; 

^amd(Z) = -^A{(Pi,<P2,-",<Pa}, (2) 

<fix t ( r j) oc exp 



<fi = 4>x t Xin : < 



\ 2 

2k 1 



Xi 




(3) 



where Xi ls the intrinsic spin function parameterized by £j, and t% is the isospin 
function which is up(proton) or down(neutron). Thus an AMD wave function 
is parameterized by a set of complex parameters Z = {X n i,£i} (n = 1 ~ 
3 and i = 1 ~ A). 

If we consider a parity eigenstate projected from a Slater determinant the 
total wave function consists of two Slater determinants, 

$(Z) = (1±P)$ AMD (Z), (4) 

where P is a parity projection operator. In case of total angular momentum 
eigenstates the wave function of a system is represented by integral of rotated 
states, 

$(Z) = Pj> IK ,$ AMD (Z) = J dnD J ^ K ,(Q)R(Q)^ AMD (Z), (5) 

The expectation values of operators by <3?(Z) are numerically calculated by a 
summation over mesh points of the Euler angles O. 

In principle the total wave function can be a superposition of independent 
AMD wave functions. We can consider a superposition of spin parity projected 
AMD wave functions P MK i&amd 7 s, 

$(Z, Z') = cP^,$amd(Z) + c'P^ ki <S>amd(Z') + ■■■■ (6) 



B. Energy variation 

We make variational calculations to find the state which minimizes the 
energy of the system; 

£ = wm (7) 

by the method of frictional cooling. Concerning with the frictional cooling 
method in AMD, the reader is referred to papers p9| , |36| . For the wave function 
3>(Z) parameterized by complex parameters Z, the time development of the 
parameters is determined by the frictional cooling equations, 

dX nk . 1 8£(Z,Z') ,.,.,1 d (*(Z)|H|*(Z)) 

-*- = (A+ *" ) a^jr = (A+ "' ) aM^ (*(z)|#(z)> ■ (n = 1 ~ 3 k = 1 ~' 4) (8) 
d£fc ,, , i ds(z,z') „ , i a ($(z)iiyi$(z)) ,, , 



with arbitrary real numbers A and [i < 0. It is easily proved that the energy 
of the system decreases with time as follows, 

d£ _^ f d£ dXi dS dXJ d£ d& d£ dg\ 

dt ~ ^ ^dx" ' ~dT + ax* ' ~dT + a&" "dt" + d£* ' ~dl) 

2n*(d£ 8£ d£ 8£\ n 

After sufficient time steps for cooling, the wave function of the minimum- 
energy state is obtained. 



C. Angular momentum projection 

Expectation values of a given tensor operator Tr(rank k) for the total- 
angular-momentum projected states are written as follows, 

(PaKk^K^Ik^) (12) 

2J 2 + 1 



8tt 2 



(j 2 M 2 kq\j 1 M 1 )J2(J2Kku\j 1 K 1 ) I "dni>^ 2 (n)($i|r*i2(n)|$ 2 ), (13) 

Ku J 



where DJ^ 1K are the well-known Wigner's D functions and R(Q) stands for 
the rotation operator with Euler angles O. In the practical calculations, the 
three-dimensional integral is evaluated numerically by taking a finite number 
of mesh points of the Euler angles 0. 

D. Simplest version of AMD for the study of nuclear structure 

In the simplest version of AMD for the study of nuclear structure, the 
ground state wave function of a system is constructed by the energy variation 
of the parity eigenstate projected from a Slater determinant. Furthermore, 
the directions of intrinsic spins of single particle wave function are fixed to be 
up and down as £j = ±5 f° r simplicity. Therefore the spin-isospin functions 
of single-particle wave function are chosen as p J, p I, n |, and n J, in the 
initial state and are fixed in the energy variation. In this case the total wave 
function of a system is parameterized only by X = {Xi, X 2 , • • • , X^} which 
are the centroids of Gaussian wave packets in the phase space, 

$(X) = (l±P)$ y4MD (X). (14) 

We regard the minimum-energy state obtained with the energy variation (de- 



scribed in II B) for the parity projected state as the intrinsic state of the 
system. In order to compare with experimental data, we project the intrinsic 
wave function to the total-angular-momentum eigenstates and calculate the 
expectation values of operators. In that sense, "the simplest version of AMD" 
stands for the variational calculations after the parity projection but variation 
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before projection (VBP) with respect to the total-angular momentum in this 
paper. In the same way as the ground state, the lowest non-normal parity 
state is calculated by energy variation for the non- normal parity projected 
state. 

E. Variation after projection 

The wave function of the system should be a total-angular-momentum 
eigenstates. We can perform energy variation after the spin-parity pro- 
jection(VAP) with the method of frictional cooling for the trial function 

^ = P J mk'^amd(Z) [Hi- 
First we make VBP calculation to prepare an initial state & am D^init) for 
the VAP calculation. We choose an appropriate K' quantum number for each 
spin parity that makes the energy expectation value for the spin parity eigen- 

State (P M K /& AMD{Z'init)\H\P M K ,& AM Di^init)) / (Pm K' i^init)\PM K'i^init)) 

minimum. K' is a component of the total angular momentum along the 
approximately principal axis on the intrinsic system. In order to obtain 
the wave function for the lowest J^ state, we perform VAP calculation for 
(P^^ A md(Z)\H\P^^ amd (Z)}/(P^^ amd (Z)\P^ k ^ A md(Z)) with 
the appropriate K' quantum number chosen for the initial state. In the VAP 
procedure, the principal z-axis of the intrinsic deformation is not assumed to 
equal with the 3- axis of Euler angle in the total angular momentum projec- 
tion. In general the principal z-axis is automatically determined in the energy 
variation. That is to say, the spin parity eigenstate P mk ,<&amd obtained by 
VAP with a given K' = ( J3) can be the state with so-called K = (J z ) quantum 
number mixing in terms of the intrinsic deformation. 



F. Higher excited states 

As mentioned above, with the VAP calculation for $(Z) = P mk ,$ A md(Z) 
of the J eigenstate with K' we obtain the wave function for the lowest J 
state, which is represented by the set of parameters Z = Z 1 .To search the 
parameters for the higher excited J^ states, the wave functions are superposed 
so as to be orthogonal to the lower states as follows. The parameters Z„ 
for the n-th J state are reached by varying the energy of the orthogonal 
component to the lower states; 

ib!7\ - pJ± dS 17\ V^ {P M K'^AMp{Z k )\P MK ,^AMP{Z)) j± j± 

n z ) ~ P M k>*amd(Z) -2^ .j± t (7 J± MP J± a j=j±-rP M Ki®AMD(Z k )• (15) 

fc=l {P M K>®AMD{Z k )\P MK ,®AMD{Z k )} 

In the present paper, we call the variational calculation after the spin 
parity projection (mentioned in previous subsection) and the calculation for 
the higher excited states described in this subsection as VAP calculations. 



G. Diagonalization in VAP 



After VAP calculations for various J^ states mentioned above, the in- 
trinsic states &\ MD , *&amd-> ' ■> ^amd> which correspond to the J^ states, 
are obtained as much as the number of the calculated levels. Finally we 
construct the improved wave functions for the J~ states by diagonalizing 
the Hamiltonian matrix (P M k ,$am D \H\P M K u^\ M D ) and the norm matrix 
(-^M K'^am d\Pm K"^am d) simultaneously with regard to (i,j) for all the in- 
trinsic states and (K',K"). In comparison with the experimental data such 
as energy levels and E2 transitions, the theoretical values are calculated with 
the final states after diagonalization. 



III. INTERACTION 



For the effective two-nucleon interaction, we adopt the Volkov No.l force 
[40] as the central force. The adopted parameters in this paper contain only 
Wigner and Majorana components. For some nuclei we have performed cal- 
culations by adding appropriate Bartlett and Heisenberg components to the 
Volkov force. However the results have proved to be not so much affected 
by the additional components except for the binding energies. Instead of the 
Volkov force, we also adopt case (1) and case (3) of MV1 force [41|, which 
contain the zero-range three-body force V^> as density dependent terms in 
addition to the two-body interaction V^ 2 ^, 



v DD = v {2) + y (3) 

yPO = (1 - m + bP a - hP T - mP^Pr) \ V A exp 



r 
ta 



+ V R exp 



r 
tr 



b = h = 0, 

V A = -83.34MeV,r A = 1.60fm, 

Vr = 104.86MeV (case I), 99.86MeV (case 3),r R = 0.82fm, 

y(3) =v ^)s(r 1 -r 2 )5(r 1 -r 3 ), 

V (V = 4000MeV fm 6 for (case 1), 5000MeV fm 6 (case 3), 

where P a and P r stand for spin and isospin exchange operators, respectively, 
and r denotes |r| = |ri — r2|. As for the two-body spin-orbit force Vls, we 
adopt the G3RS force @ ; 

V LS = {ui exp (- K/ r 2 ) + u n exp (-K H r 2 ) } P( s O)L • (Si + S 2 ), 
.d 



L = r x 



di 



ki = 5.0fm" 2 ,K// = 2.778fm -2 , 



with P( 3 0) denoting the projection operator onto the triplet odd( 3 0) two- 
nucleon state. The Coulomb interaction is approximated by a sum of seven 
Gaussians. 



(16) 

(17) 



(18) 
(19) 



IV. STUDY WITH SIMPLEST VERSION OF AMD 

Since the wave function should be a total-angular-momentum eigenstate, 
it is expected that the VAP calculation gives better results than the sim- 
plest version of AMD. It is, however, not easy to perform VAP calculation 
because of the three-dimensional integral for the total-spin projection which 
is evaluated by taking a large number of mesh points of the Euler angles. In 
order to study systematically the structures of ground states of light nuclei 
covering from the ordinary region to the exotic unstable region, we perform 
the simplest version of AMD calculations (energy variation for the parity pro- 
jected state) for even-odd, odd-even and even-even isotopes of Li, Be, B and 
C. The obtained states are projected to the total- angular-momentum eigen- 
states in calculating the expectation values to compare the results with the 
experimental data. Fortunately, in many nuclei with some exceptions, the 
obtained structures are not so much different from the ones obtained with 
VAP calculations. The results with VAP will be shown in the next section. 

A. Results 

In this section the theoretical results of the simplest version of AMD are 
compared with the experimental data. 

We have used the Volkov force with Majorana parameter m = 0.56 ((a)), 
and the case 1 of MV1 force with m = 0.576 ((b)) for the Li and Be isotopes. 
We have also adopted the MV1 force with mass dependent m parameters ((c)) 
for B isotopes, and the MV1 force with m = 0.63 ((d)). The adopted param- 
eters are listed in Table | and |TJ. The details are explained in each place. The 
optimum width parameters v are shown in Table |TO| , |IV|,|V|. The expectation 
values for the operators of observable quantities are calculated by projecting 
the intrinsic states obtained with the simplest version of AMD into the eigen- 
states of parity and total-angular-momentum. Most theoretical values except 
for the energy levels are obtained without the mixing of if -quantum which 
is the component of the spin J along the approximate principal axis of the 
intrinsic system. Instead, we choose an appropriate K which gives the mini- 
mum energy for each spin of a system. We have diagonalized the Hamiltonian 
matrix with respect to the K quantum number within the spin J projected 
states to calculate energy levels, and found only small mixing of the K quan- 
tum number which implies that K is approximately a good quantum number 
in the lowest J states projected from the intrinsic states. 

1. Energies 

Figure [I] shows the binding energies of the ground states of Li and Be 
isotopes. With both of the interactions (a) and (b), the binding energies of Li 
and Be isotopes are qualitatively reproduced. The calculated result of n Be is 
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TABLE I. The adopted interaction parameters. 



(b) 
(d) 
(e) 



Volkov No.l 
MV1 case(l) 
MV1 case(l) 
MV1 case(l) 



m = 0.56, b = h 

m = 0.576, b = h 

m = 0.63, b = h 

m = 0.336, b - 







0.2, h = 0.4 



UI = -u n = 900 MeV 
UI = -u n = 900 MeV 
UI = -u n = 900 MeV 
uj = -u H = 1500 MeV 



TABLE II. The mass dependent interaction parameters of the force (c) adopted for B isotopes . 



n B, 13 B 


MV1 case(l) 


m = 0.576, b = h = 


15 B 


MV1 case(l) 


m = 0.63, b = h = 


17 B, 19 B 


MV1 case(l) 


m = 0.65, b = h = 



UI = -u u = 900 MeV 
UI = - UII = 900 MeV 
UI = - Uu = 900 MeV 



TABLE III. The width parameters v of Gaussian wave packets adopted in the calculations of 
the simplest version of AMD for Li and Be isotopes. The optimum width parameter is chosen for 
each nucleus and each set of interactions. 





width parameter 


v (fm 2 ) 


interaction 


(a) 


(b) 


6 Be(+) 


0.215 


0.195 


7 Be(-) 


0.230 


0.200 


8 Be(+) 


0.250 


0.205 


9 Be(-) 


0.245 


0.195 


9 Be(+) 


0.235 


0.200 


10 Be(+) 


0.230 


0.190 


10 Be(-) 


0.225 


0.190 


n Be(-) 


0.220 


0.180 


n Be(+) 


0.220 


0.180 


12 Be(+) 


0.215 


0.175 


12 Be(-) 


0.210 


0.180 


13 Be(+) 


0.205 


0.170 


14 Be(+) 


0.210 


0.170 


7 Li(") 


0.230 


0.200 


9 Li(+) 


0.210 


0.180 


H Li(+) 


0.195 


0.170 
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TABLE IV. The optimum width parameters v of Gaussian wave packets adopted for B isotopes. 









Width parameter 


v (fm 2 ) 






Interaction 


n B 


13 B 


15 B 




17 B 


19 B 


m=0.576 


0.185 


0.175 


0.175 




0.165 


0.160 


777=0.63 


0.170 


0.160 


0.155 




0.150 


0.145 


777=0.65 


— 


0.155 


0.150 




0.150 


0.135 



TABLE V. The optimum width parameters v of Gaussian wave packets adopted for C isotopes. 







width parameter v (fm 2 ) 


interaction 


(b) 


(d) 


9 c(-) 


0.18 


0.170 


10 C(+) 


0.19 


0.180 


n c(-) 


0.19 


0.175 


12 C(+) 


0.19 


0.175 


13 c(-) 


0.18 


0.170 


14 C(+) 


0.18 


0.165 


15 C(+) 


0.17 


0.160 


16 C(+) 


0.17 


0.160 


17 c(+) 


0.17 


0.155 


18 C(+) 


0.17 


0.150 


19 C(+) 


0.16 


0.150 


20 C(+) 


0.16 


0.145 


22 C(+) 


0.16 


0.140 
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FIG. 1. The binding energies of Li and Be isotopes. The circles(triangles) are the theoretical 
values with the interaction (a) ((b)). The square points indicate the experimental data. 



> 

CD 



1— 

CD 
c 
LU 

c 



CO 




7 9 11 

mass number A 



8 10 12 

mass number A 



the binding energy of the lowest \ state, though the normal-parity \ state 
is lower than i state in these calculations. Detailed discussion of the energy 
levels and the parity of the ground state of n Be will be given later. 

The binding energies of B and C isotopes are presented in Fig.^ and Fig.||, 
respectively. For B isotopes, the figure shows theoretical results with the MV1 
force with m = 0.576, m = 0.63, and m = 0.65. In Fig. y, results of C isotopes 
with Volkov(m = 0.60) and MV1 force with m = 0.576 and m = 0.63 are 
shown. For 15 C the energies of 5/2 + states are shown, although the ground 
state should be l/2 + . In both B and C isotopes, MV1 force with m = 0.576 
reproduces well the experimental data. 

The energy levels of Li, Be, B, and C isotopes are displayed in Fig.[|, 
Fig.||, Fig.[| and Fig.|?]. The adopted interactions are explained in the figure 
captions. The second J states are obtained by diagonalizing Hamiltonian 
with respect to K quantum numbers. In many nuclei, the theoretical values of 
rotational bands such as excitation energies and spin sequence in low energy 
region well correspond to the experimental data. It means that many low- 
lying levels are approximately described as the rotated states of the intrinsic 
states obtained by the simplest AMD calculations which are VAP for parity 
projection but VBP for total-spin projection in the framework of AMD. Com- 
paring the level spacing calculated with and without the three-body force, it 
is found that the states obtained with the three-body force have smaller level 
spacing, that is, larger moment of inertia. In many cases we found that the 
VBP calculation gives smaller level spacing than the VAP calculation. 

The energy levels depend on the parameter m of the Majorana exchange 
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FIG. 2. The binding energies of B isotopes. The circles are the theoretical values calculated 
with the MV1 force with m = 0.576, 0.63, and 0.65. The square points indicate the experimental 
data. 
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term and the presence of the three-body force. In general the results without 
the three-body force overestimate the excitation energies of non-normal par- 
ity states in most nuclei. The differences of the excitation energies between 
the normal parity and the non-normal parity states are improved in the re- 
sults with the three-body force. In general the non-normal parity states have 
wider extension of the density distribution and therefore they feel relatively 
weaker repulsive force due to the three-body terms. It is the reason why the 
calculations with the three-body force give smaller excitation energies for the 
non-normal parity states. Even in the results with the three-body force, how- 
ever, the non-normal parity l/2 + state is still higher than the l/2~ states in 
n Be which has been known to have abnormal parity of the ground state l/2 + . 
Careful choice of interactions and the improvement of the wave function is 
necessary to reproduce this feature of the parity inversion. For instance, VAP 
calculations with appropriate interaction parameters succeed to obtain the 
lower energy of l/2 + state than the one of l/2~ state. The detail will be 
mentioned later. 

The energy levels of odd-even and even-odd nuclei are also sensitive to 
the strength of the spin-orbit force. The other set of interaction parameters 
(f ) explained in the figure caption of FigJEI give rather good results for energy 
levels of n B. 

In 13 C and 15 C the lowest positive parity states are known to be l/2 + 
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FIG. 3. The binding energies of C isotopes. The lines are the theoretical values calculated 
with the MV1 force with m = 0.576 and 0.63. The figure also shows the calculations with the 
Volkov No.l force with m = 0.6. The square points indicate the experimental data. 
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states. However, the calculated intrinsic states of positive parity with 
MVl(m = 0.576) and Volkov force(m. = 0.60) contain little component of 
l/2 + states. Adopting other interaction parameter (e) in Table Q: MV1 
force with m = 0.336, Bartlett and Heisenberg components b = —0.2, 
h = 0.4 and the slightly stronger spin-orbit force with the magnitude of 
ui = —Uii = 1500MeV, we obtained the l/2 + component whose energy is 
still rather higher than the one of 5/2 + state. In the result with every inter- 
action, the main component of the intrinsic state obtained with VBP calcu- 
lation is 5/2 + , therefore, the variation affects to minimize the energy of 5/2 + 
state. For the lowest l/2 + states of 13 C and 15 C, variations after spin-parity 
projection are useful instead of VBP calculations. 

2. Radii 

Figure || shows the radii of Li and Be isotopes. Dashed lines are the theo- 
retical root-mean-square radii by AMD with the three-body force and dotted 
lines show results without three-body force. Square points indicate the ex- 
perimental radii reduced from the experimental data of the interaction cross 
section. The AMD calculations with the three-body force seem to qualita- 
tively agree with the observed radii except for the very neutron-rich nuclei. 
The theory underestimates the extremely large radii of 11 Li, n Be and 14 Be 



14 



FIG. 4. The energy levels of Li and Be isotopes. The calculations with the interaction (a) 
Volkov No.l force (m = 0.56) and (b) MV1 force (m = 0.576) are compared with the experimental 
data. 
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FIG. 5. The energy levels of B isotopes. The adopted interactions are the case (b) MV1 
force (m = 0.576). The results for n B with the MV1 (m = 0.56) and the spin-orbit force 
(uj = — u//=1500MeV) are also shown (f). 
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which are considered to have the neutron halo structures. For reproduction 
of such large radii due to the halo structures, improvement of the wave func- 
tion and careful choice of interactions should be important. The radius of 
the positive parity state of 12 Be calculated with the present simple AMD cal- 
culation is smaller than the experimental datum because the obtained state 
has the closed neutron-p-shell structure. However VAP calculation with a set 
of interaction parameters which reproduces the parity inversion of the n Be 
ground state gives a ground state of 12 Be as 2Tioj state with 2 particles in sd 
shell and 2 holes in p shell in neutron configuration, whose radius is as large 
as the experimental data. 

In Fig. 0, theoretical radii of B isotopes are compared with the experi- 
mental interaction radii. The triangles connected with the dotted line indicate 
the AMD results by the use of the fixed Majorana parameter m = 0.576 (in- 
teraction(b)). The solid line shows the results by the interaction (c) with a 
mass-dependent Majorana parameter; m =0.576, 0.576, 0.63, 0.65, and 0.65 
for n B, 13 B, 15 B, 17 B, and 19 B, respectively. Considering that the use of 
larger m value for the heavier nuclear system is generally reasonable, it is not 
unnatural to adopt the mass-dependent m values adopted here. The results 
with mass-dependent Majorana parameter reasonably fit to the experimental 
data. 

The radii of C isotopes are presented in Fig.|l0| The theoretical results 
are calculated by the use of the MV1 force (the solid line for m = 0.576 and 
the dotted line for m = 0.63). As seen in the figure, the recent experimental 



data of the interaction radii of C isotopes [43] are found to be consistent 
with our theoretical predictions [H. The radii of C isotopes have a kink 



at C and increase as the neutron number becomes larger in the neutron- 
rich region N > 8. It is easy to quantitatively fit the theoretical values to 
the experimental ones except for the valley at n C by using mass-dependent 
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FIG. 6. The energy levels of C isotopes. The theoretical values are obtained with the Volkov 
force (m = 0.6) (f) and with the interaction (b) MV1 force (m = 0.576). For the positive parity 
states of 13 C and 15 C, the results with the interaction (e) are presented instead of the ones with 
the interaction (b). 
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FIG. 7. The energy levels of C isotopes. (See the figure caption in Fig.| 
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FIG. 8. The root-mean-square radii of Li and Be isotopes. The solid lines and the dotted 
lines are the AMD results calculated with the interaction (a) Volkov force with m = 0.56, and (b) 
MV1 force with m = 0.576. The square points indicate the experimental data deduced from the 
interaction cross sections |5|. 
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FIG. 9. The root-mean-square radii of B isotopes. The triangles connected with the dotted line 
indicate the AMD results by the use of interaction (b) the fixed Majorana parameter m = 0.576. 
The solid line shows the results by the interaction (c) with a mass-dependent Majorana parameter. 
The square points indicate the experimental data deduced from the interaction cross sections fq]. 
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FIG. 10. The root-mean-square matter radii of C isotopes. The solid and dashed lines are the 
theoretical results by the use of interaction (b) m = 0.576 and (d) m = 0.63. The experimental 
data deduced from the interaction cross sections are displayed with the square points [43|. 
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Majorana parameter in a similar way to the case of B isotopes. In AMD 
results, we did not find any reason for the small radius of n C. A detailed 
discussion of the radii of the neutron-rich C isotopes will be given later in 



subsection IV B 4 about the neutron-skin structure. 



3. Magnetic moments 

Magnetic dipole moments by the of AMD wave function are obtained by 
numerically calculating the expectation values of the magnetic dipole operator 
p, z by the spin-parity projected states with the highest z-component of the 
spin, 

where jX z is the z component of the operator of the magnetic dipole moment 
fi = g sn Sn + 9spSp + ddp with the bare ^-factor g sp = 5.58, g\ = 1.0 for protons 
and g sn = —3.82 for neutrons. We choose an appropriate K quantum number 
which minimizes the energy of the state P JK &. Figure 11 shows the magnetic 
dipole moments fi of odd-even Li isotopes, even-odd Be isotopes, and odd- 
even B isotopes. The theoretical results of AMD calculations agree with the 
experimental data for many nuclei very well. It should be emphasized that 
the AMD method is the first framework which has succeeded in reproducing 
the magnetic dipole moments for these isotopes systematically. 

The dependence of the \x moment on the neutron number in Li isotopes 
and in B isotopes is closely related to the structure change with the increase 
of the neutron number. Therefore these data of \x moments should carry 
important informations about the nuclear structures. We will give detailed 
discussions about the correlation of the nuclear structure and the observed 
electromagnetic properties in later section. 

The theoretical value of the /x moment of 11 Be(l/2 + ) is sensitive to the 
strength of the spin-orbit force. In the case of MV1 force with m = 0.576, 
the fj, moment of 11 Be(l/2 + ) calculated with the strength of the spin-orbit 
force ui = uji = 900MeV is —1.9 /j,n which is as much as the Schmidt value, 
while with the stronger spin-orbit force ui = uu = 1200MeV the /i moment 
is —1.2 //TV- In the case of l/2 + state of 11 Be, the \i moment directly reflects 
the spin configuration of the last 7-th neutron. Since the strength of the spin- 
orbit force affects the orbit of the last neutron in sd shell and also the core 
excited component of the neutron p 3 / 2 closed shell, it is natural that the the 
H moment of 11 Be depends on the strength of the spin-orbit force. 

On the other hand, in the case of /i moments of odd-even Li and B isotopes 
the calculated results do not depend on the interaction parameters so much as 
the one of 11 Be because the main contribution to [i moments originates from 
the the spin configuration of the valence proton in the P3/2 orbit. A slight 
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(20) 



FIG. 11. Magnetic dipole moments fj, of even-odd Be and odd-even Li and B isotopes. The 
calculations with the interaction (b) are compared with the experimental data. 
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dependence of \x moments of 9 Li and the mirror nucleus 9 C on the spin-orbit 
force has been discussed in Ref . [ 15 1 . 



In Table |V| the \i moments of even-odd C isotopes are presented. The 
theoretical results are calculated with the interaction (b) MV1 force with m = 
0.576 and the spin-orbit force with ui = —uu = 900MeV except for positive 
parity states of 13,15 C. The positive parity states of 13 > 15 C are calculated with 
the interaction (e) MV1 force with m = 0.336, b = —0.2 and h = 0.4 and the 
spin-orbit force with the magnitude of uj = —ujj = 1500MeV which is the 



same interaction as mentioned in IV A 1 . The theoretical results well agree to 
the experimental data. 



Electric quadrupole moments and B(E2) 



It should be pointed out that we can describe electric properties such as 
quadrupole moments and E2 transition strength B(E2) by using not effective 
charges but the bare charges for protons and neutrons in the AMD framework. 
It is because the drastic changes of proton and neutron structures are directly 
treated in the framework. Here, we just present the theoretical results of 
electric quadrupole moments and B{E2) comparing with the experimental 
data. We will give detailed discussions on the relation between observable 
electric properties and the intrinsic structures later in IV B| . 

Although the intrinsic structures of many nuclei are qualitatively not so 
much sensitive to the adopted interaction, the calculated Q moments depend 
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FIG. 12. Electric quadrupole moments Q of Li, Be, and B isotopes. The theoretical values 
and the experimental data for the ground 3/2~ states of odd-even Li and B isotopes, and 9 Be 
are shown. Solid lines (A) and triangle points indicate the calculated results using the MV1 force 
with m = 0.576 (interaction (b)) for Li and Be, and using the mass-dependent m parameters 
(interaction (c)) for B isotopes. The dashed line (B) indicates the results of B isotopes with 
Bartlett and Heisenberg terms described in the text. The point (C) is the calculation of 'Li with 
the wave function improved by taking account for the long tail of the relative motion between 
clusters a and t. The experimental data are quoted from Refs. flS,lC,46| 
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TABLE VI. Magnetic dipole moments and electric quadrupole moments of C isotopes. We 
adopt the interaction (b) except for positive parity states of 13 C and 15 C. The interaction (e) is 
used for positive parity states of 13 C and 15 C. The experimental data are quoted from Ref. fl2jl a , 
and Ref. p| b . 
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on the adopted interaction with or without a density-dependent term and also 
depend on the Majorana parameter m, because the radii are sensitive to these 



interaction parameters. Figure |12| shows the electric quadrupole moments of 
Li, Be, and B isotopes. The solid lines and the triangle points (A) indicate 
the AMD calculation of Li and Be isotopes by the use of MV1 force with m = 
0.576 (interaction (b)), and that of B isotopes with the mass-dependent m 
parameters (interaction (c)) which reproduce the observed radii of B isotopes. 
The calculations agree to the experimental Q moments systematically . 

In order to check the interaction dependence of Q moments we also show 
the theoretical values for B isotopes calculated by adopting the interaction 
with additional Bartlett and Heisenberg terms as b = —0.2 and h = 0.4 in 
Fig.^J (dashed line). Here the Majorana parameters are changed from m into 
m = M + -r(b — h) so as to give as the same binding energies for a nuclei as the 
ones obtained with the interaction (c) of no Bartlett and Heisenberg terms. 
Comparing the solid line and the dashed line for B isotopes, it is found that 
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the additional components do not give significant effects on Q moments. 

Calculations with the interaction (b) fit well to the experimental Q mo- 
ments of Li and Be isotopes except for 7 Li. The Q moment of 7 Li is underesti- 
mated by theory. By improving the wave function of 7 Li in the following way, 
we have obtained the large Q moment of Li indicated by a point (C) in Fig. 



12 which is as much as the experimental data. As we explain the calculated 
intrinsic structure later in detail, the AMD wave function of 7 Li has proved 
to have the well-developed cluster structure of a + t. When the clustering 
is well developed, the relative wave function between clusters spreads out to- 
ward the outer spatial region resulting in a long tail. However, since the single 
nucleon wave function of AMD is a Gaussian wave packet, the relative wave 
function between clusters is also necessarily a Gaussian wave packet. Because 
of the lack of the outer tail part of the relative wave function between a and 
t , the quadrupole moment of 7 Li is underestimated by the simplest AMD 
wave function. Therefore we have improved the inter-cluster relative wave 
function of the AMD by superposing several AMD wave functions which are 
written asa + t clustering states with different distances between the centers 
of two clusters. Superposition of the spin-parity eigen states projected from 
these wave functions has been made by diagonalizing the total Hamiltonian. 
The improved wave function has proved to reproduce the electric quadrupole 
moment well as seen in Fig.[TT^, 

Also in the case of C, the calculations with the interaction (b) reasonably 
agree with the observed Q moments of n C and 12 C. As seen in the radii 
of unstable nuclei with neutron halo and also in the Q moments of 7 Li, the 
simple AMD wave function is not sufficient to describe the outer tail of wave 
function because of the Gaussian form. If the proton-rich nuclei such as 9 C 
have the outer tail of the valence protons, the theoretical predictions may 
underestimate the Q moments of these nuclei. 



Figure 13 shows the E2 transition strength. The theoretical values are 
calculated by using the interaction (b), MV1 force with m = 0.576, except 
for the data of the positive parity states of 15 C. Theoretical values well agree 
with the experimental data. B(E2;l/2~ — ► 3/2 _ ) in 7 Li is underestimated 
by simple AMD wave function because of the lack of outer tail of the inter- 
cluster (a-t) relative wave function. The strength can be reproduced by the 
theoretical results with the improved wave function described above. In the 
simple AMD results the strength B(E2; 2+ — ► + ) of 12 Be is calculated to be 
smaller than the one of 10 Be. However a larger strength of B(E2; 2 + — > + ) 
of 12 Be than the one of 10 Be is predicted by the VAP calculations with a set 
of interactions which gives a largely deformed ground state with 2p — 2h in 
the neutron orbits. 
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FIG. 13. El transition strength in Li, Be, B, and C isotopes. Square points indicate the 
experimental data and triangle points are the AMD calculations with the interaction (b) MV1 
m = 0.576. The initial and final states are the lowest J 71 " levels except for the levels with subscripts 
in n B. The cross point for B(E2 : 1/2 — ► 3/2) of 7 Li is the calculation with the wave function 
improved by taking into account the long tail of the relative motion between a and t clusters. 
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FIG. 14. The density distributions of the normal-parity states of Li isotopes. The density of 
the intrinsic wave function before the parity projection is integrated along the line parallel to an 
appropriate principal axis as explained in the text. The AMD wave functions are obtained with 
the interaction (b). 
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B. Discussion 

In this section we make systematic study of the intrinsic structures. Af- 
ter describing the deformations and clustering structures of Li, Be, B and 
C isotopes, we discuss the effect of the intrinsic structure on the observable 
quantities to deduce the informations of the nuclear structure from the ex- 
perimental data. We notice some interesting features such as the effect of 
clustering structures on the electromagnetic properties, the opposite defor- 
mation between protons and neutrons, and the neutron skin structures. 

1. Shapes and clustering structure 



By analyzing the intrinsic wave function, we discuss the shapes, defor- 
mations and clustering aspects of Li, Be, B and C isotopes. First we show 
the density distributions of the intrinsic states of Li, Be, B and C isotopes 
in the Figs. 14, 15, 16, [IR and IS. In drawing the figures, the density 



of each intrinsic wave function before parity projection is projected onto an 
X-Y plane by integrating out along the line parallel to the Z axis. Here 
X, Y, Z axes are chosen so as to be (X^ x ?) > ^liVi) > (J2i z i) an d 
{J2i x iVi) = (J2illi z i) = (J2i z i x i) = 0- We see the drastic structure change 
along the increase of the neutron number in each series of isotopes. 



In the results of Li isotopes (Fig. 14), it is easily seen that the 7 Li system 
has the largest deformation with the a+t clustering structure. 9 Li also has a 
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FIG. 15. The density distributions of the normal-parity states of Be isotopes. The adopted 
interaction is the case (b). 
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deformed shape, though the deformation is not as large as the one in 7 Li. The 
ground state of n Li has an almost spherical structure that can be expressed 
by a shell model wave function with the closed neutron p shell. 

As for the Be isotopes (Fig. 15 and[l^), separated two pairs of protons are 
found in the proton density. It means that 2a core exists in all the heavier Be 
nuclei than 7 Be. The development of clustering in Be isotopes is estimated 
by the relative distance between two pairs of protons, which we show in Fig. 



17. In the non- normal-parity states we find that the la clustering is largest 
in 9 Be as already well known. In the normal-parity states of Be isotopes, the 
clustering becomes weaker and weaker with the increase of neutron number up 
to 12 Be with neutron magic number N = 8. In 14 Be the clustering develops 
again and the deformation becomes larger than 12 Be. In the non-normal- 
parity states of Be isotopes, there are many exotic structures with developed 
clustering structures and larger deformations than the normal parity states 
(Figjlq). These largely deformed states give rise to the rotational bands which 
well agree to the experimental data of energy levels of the non-normal-parity 
states in 9 Be and 10 Be. The ground state of n Be is known to be a non- 
normal state with l/2 + . The calculated positive-parity state of 11 Be which 
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FIG. 16. The density distributions of the non-normal-parity states of Be isotopes. The adopted 
interaction is the case (b). 
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FIG. 17. The distance R pp between two pairs of protons in the intrinsic states of the normal 
parity states (solid line) and the non-normal parity states (dotted line) of Be isotopes. In Be nuclei 
with JV > 4 the distance R pp is considered to be the inter-cluster distance between 2a. 
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corresponds to the ground l/2 + state has the developed prolate deformation 
as large as the normal-parity state of 9 Be (see Fig. |l7]). The deformation in 
the positive parity state of n Be is considered to be one of the reasons for the 
parity inversion of the ground state. 

Here we stress again the possibility of the abnormally deformed ground 
state of 12 Be. With the simple AMD calculations, the obtained normal-parity 
state of 12 Be has the closed p-shell structure with a spherical shape of neu- 
trons. However, the VAP calculations with the set of interaction which re- 
produces the abnormal spin parity l/2 + of the ground state of n Be suggest 
that the ground state of 12 Be is a 2hu> state with 2 neutrons in sd shell in the 
language of a simple shell-model. In that case, the ground state of 12 Be has a 
large prolate deformation with a developed clustering structure, and instead, 
the spherical p-shell closed-shell state is found in the second 0^ state. In 
the later section on VAP, we will mention the details about the structure of 
excited states of neutron-rich Be isotopes. Although the lhco and 2Tuj states 
are suggested to be the ground states of Be isotopes in N > 7 region, in this 
section based on the simple AMD results we discuss the so-called OHlo states 
of normal-parity states and Vhuj states of the non-normal parity states which 
are expected to be the ground or low excited states. 

Also in the AMD results of B isotopes shown in Fig. 



18, the drastic 



structure change with the increase of the neutron number is found. The 
total matter density p indicates the deformed state with a three-center-like 
clustering in n B, while the nucleus 13 B with a neutron magic number N = 8 
has the most spherical shape among B isotopes. It is very interesting that 
in the neutron-richer region, 15 B, 17 B and 19 B, the clustering structure with 
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FIG. 18. The density distributions of the normal-parity states of B isotopes. The intrinsic 
wave functions are calculated by the use of the interaction (c) with the mass-dependent Majorana 
parameter. 
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FIG. 19. The density distributions of the normal-parity states of C isotopes. The intrinsic 
wave functions are calculated with the interaction (b) except for 15 C. For 15 C, we obtain the 
intrinsic state with the interaction (e). 
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prolate deformation develops again. In the right column for neutron density of 
Fig. 18 we can see the neutron structures. In 11 B, six neutrons have an oblate- 



deformed distribution, while eight neutrons in 13 B constitute the closed shell. 
On the other hands, ten neutrons in 15 B posses a large prolate deformation. 
It is consistent with the familiar features that the ordinary nucleus 12 C with 
N = 6 has an oblate shape , the nucleus 16 with N = 8 is spherical and 
20 Ne with N = 10 posses a large prolate deformation. The neutron densities 
in 17 B and 19 B are found to have largely prolately deformed structures. The 
prolate deformation in the system with N = 14 is not an obvious feature, but 
is a characteristic seen in neutron-rich B isotopes in which five protons prefer 
a prolate shape. Generally speaking proton structure follows the change of 
neutron deformation in B isotopes. It means that in the neutron-rich region 
the proton density with two clusters stretches outward as the neutron number 
goes up toward the neutron-drip line. As a result we find that the two-center- 
like clustering develops more largely in 17 B and most largely in 19 B. It is 
interesting that the theory predicts the development of the clustering in B 
isotopes near the neutron-drip-line as well as Be isotopes. The present results 
for the development of clustering are consistent with the previous work by 
Seya et al. |47| where they calculated B isotopes by assuming the existence 
of a 2a core. It should be pointed out that the present work by AMD is 
the first calculation which predicts the clustering structure in neutron-rich 
Be and B isotopes without a priori assumption of the cluster cores as far as 
we know. We consider that these clustering structures in the unstable nuclei 
near the neutron drip-line consist of cluster cores and surrounding neutron 
cloud because the valence neutrons are weakly bound in these nuclei. In that 
sense the clustering structures seen in neutron-rich region should be different 
from the well-known clustering structure in ordinary nuclei without valence 
nucleons around clusters. 

In contrast to Be and B isotopes, clustering structures are not found in 
the normal-parity states of neutron-rich C isotopes. The results show the 
general tendency of the oblately deformed proton density in the C isotopes. 
According to the AMD calculations, the well-known 3a clustering in 12 C has 
been checked with the framework free from the assumption of the existence of 
clusters. The figure for neutron density (right column in Fig.|i~9|) presents the 
drastic change of neutron structure with the increase of the neutron number. 
In the neutron-rich region the neutron density stretches widely in outer region. 
However proton density does not change so much in spite of the drastic change 
of neutron structure and remains in the inner compact region. The stability of 
the proton structure is a characteristic of the neutron-rich C with six protons. 
As a result the neutron skin structure may appear in the neutron-rich C. We 
will make more quantitative discussions of the neutron skin in C later. In 
the non-normal parity states of the proton-rich C isotopes, there are many 
exotic shapes with the large deformation. The well-developed 3a clustering 
of 12 C seen in the negative parity states constructs a rotational K = 3~ band 
which well corresponds to the lowest negative-parity band observed in the 
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experimental levels. 

For the sake of the systematic study of deformations of proton density and 
neutron density, here we explain the deformation parameters defined by the 
moments (x 2 ), (y 2 ), and (z 2 ) for the intrinsic AMD wave functions as follows; 
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Here the directions x, y, and z are chosen so as to satisfy the relation 
(x 2 ) < (y 2 ) < (z 2 ). The deformation parameters f3 7 calculated for pro- 
tons and for neutrons are displayed in Fig. 20. The Figures are for the 



normal-parity states of Li, Be, B and C isotopes with the even neutron num- 
ber obtained with the interaction (b) m = 0.576 for the nuclei with N < 8 
and m = 0.63 for the nuclei with N > 8. The behavior of the deformation 
parameters is not so sensitive to the Majorana parameter. It is found that 
the neutron shape changes rapidly with the increase of neutron number in 
all the series of isotopes. In the region N < 12, the main feature of defor- 
mations of the neutrons, prolate or oblate shape, is dominated only by the 
neutron number. It means that the neutron shape is not sensitive to the 
proton number in the light region. The neutron density in the system with 
N = 2 has a spherical shape. In N = 4 system the neutron density deforms 
prolately, and in the case of N = 6 the neutrons prefer oblate deformation. In 
N = 8 the neutron shape becomes almost spherical due to the closed neutron 
p shell. When the neutron number becomes 10, the prolate neutron deforma- 
tion appears again. As for the proton deformation, it depends on the proton 
number. In the system with Z = 2, Be isotopes, the proton density prefers 
the prolate deformation. The degree of prolate deformation of Be isotopes 
changes following the drastic change of the neutron deformation. Especially 
the prolate deformation of protons is well developed in the system with the 
prolately deformed neutron density. Contrary to Be isotopes, in the case of C 
isotope with Z = 6, the proton density always prefers the oblate deformation. 
The deformation parameter for proton is stable in spite of the change between 
prolate and oblate shapes of the neutron density. In the case of Z = 3 and 
Z = 5, Li and B isotopes, the proton shape depends on the neutron number 
so as to follow the deformed mean field given by the neutron density. In the 
system with the heavier neutron number such as N = 14, the neutron shape 
possesses both characters, which are seen in the oblate neutron shape in 20 C 
and in the prolate deformation in 19 B. In such a case with a few choice of the 
neutron shape, the neutron shape is determined by the proton number. 

As mentioned above, in the very light region the neutron deformation is 
dominated by the neutron number, and the proton shape is basically deter- 



(21) 
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FIG. 20. The deformation parameters of the intrinsic states of odd-even Li and B isotopes 
and even-even Be and C isotopes. The definitions (3 and 7 are described in the text. The adopted 
interactions are the force (b) MV1 with m = 0.576 for the nuclei N < 8 and the MV1 force with 
m = 0.63 for the heavier nuclei N > 8. 
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mined by the proton number. The dependence of neutron(proton) shape on 
the neutron (proton) number is consistent with the ordinary understanding 
of the shell effect for the stable nuclei. That is to say the spherical shape 
is seen when the neutron number equals the magic number M, the prolate 
deformation at M + 2, and oblate one at M — 2. However one of the new 
features found in this study of the light unstable nuclei is that the proton 
shape and the neutron shape do not necessarily correlate together in light 
region. As a result interesting phenomena such as the opposite deformation 
between protons and neutrons may occur in unstable nuclei. For instance, 
C nuclei prefer oblate deformation of protons while neutrons tend to deform 
prolately when the neutron number equals to 4 and 10. Therefore 10 C and 
16 C may have the opposite deformation between protons and neutrons. The 
detail of this problems in proton-rich C will be discussed later. The other 
interesting feature is the large deformation in the nuclei with prolate protons 
and prolate neutrons. The developments of deformations in 8 Be, 14 Be, 15 B, 
17 B, and 19 Be are understood as follows. Once a prolate shape of protons is 
chosen by the proton number, the prolate deformation is enhanced and the 
clustering of protons is developed by the neutron deformations if the neutrons 
deform prolately. 

In order to analyze clustering development quantitatively, we calculate the 
expectation values for the total number of the oscillator quanta. In general 
if the clustering structure develops the wave function of the system contains 
the components of a large amount of the orbits in higher shells, therefore 
the the expectation values of the oscillator quanta become larger. On the 
other hands a small value of the oscillator quantum indicates that the state 
is almost written by the shell-model states in the Ohuj configurations but the 
spatial clustering does not exist. We introduce the value AN p (AN n ) which 
stands for the deviation of the proton (neutron) orbits in the AMD wave 
function from those in the Ohcv states of the harmonic oscillator shell model 
basis, 

AM — ( P Jlk®AMD I N° p I P^^amd) N min (nn\ 

AAT— (P'mk^amd I N° p I P mk <$>amd) N min ( O o\ 

(24) 

where N° p and N° p are the oscillator quantum number operators and N p nm 
(N™ m ) are the minimum values of oscillator quantum numbers for protons 
(neutrons) given by the OHuj state. We choose the same width parameters 
for N° p and N° p as the width of the Gaussians in the AMD wave functions 
for simplicity. In Fig. ^l], AN P and AN n for the normal parity states are 
displayed as a function of neutron number. We can estimate the clustering 
development by the deviation of proton orbits AN p from the Ohuj state. As 
is expected, the dependence of AN p of Be isotopes on the neutron number is 
found to be very similar to the that of the inter-cluster distance between 2a 
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FIG. 21. The deviation of the proton orbits and neutron orbits from the harmonic oscillator 
shell- model wave functions in the normal-parity states obtained with AMD calculations. The 
expectation values of the total number of the oscillator quanta are calculated with the spin-parity 
eigen states projected from the intrinsic AMD wave functions. The detail is mentioned in the text. 
The interaction (b) MV1 with m = 0.576 is used in the region N < 8, and the MV1 force with 
m = 0.63 is adopted in the region N > 8. 
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which has been already shown in Fig. 17. In the Li, Be, and B isotopes the 



neutron number dependence of AN p is qualitatively similar to the behavior of 
A N n for neutrons. It is because the development of the clustering structure in 
these isotopes is sensitive to the neutron structure determined by the neutron 
number. In the Be and B isotopes, the increase of AN p in the region N > 8 
indicates that the clustering develops as the neutron number increases toward 
the neutron-drip line. In all the nuclei with N = 8, the AN n and also AN P 
are almost zero, which stands for that the states are approximately same with 
the shell-model states and can be written by the major shell orbits. We stress 
again that we analyzed the structure change of the normal-parity states where 
the main components are Ohio states even if they are not necessarily ground 
states. The effect of the neutron magic number N = 8 is clearly seen in the 
valley at N = 8 in the variation of AN P . In the case of the C isotopes, the 
large AN p at N = 6 indicates the developed clustering structure in 12 C, while 
the small values of AN p in the region N > 8 are because of the disappearance 
of clustering in the neutron-rich C as already mentioned. 



2. Effects of intrinsic structures on electromagnetic properties 

Here we consider the effect of clustering structures on the observable elec- 
tric and magnetic moments by analyzing /j, and Q moments of the ground 
states in Li and B isotopes. As shown above the simple AMD calculations 
well agree to the experimental data of the electromagnetic properties such as 
the magnetic dipole moments \x and the electric quadrupole moments Q, the 
strength of transitions. Roughly speaking, the electromagnetic properties of 
the odd-even nucleus reflect the orbit of the last valence proton. In that sense 
the last proton in p-shell may dominate the moments in the odd-even Li and 
B isotopes, however, jjl and Q moments in these nuclei shift as functions of 



the neutron number (Fig. 11 and 12). The N dependence of the experimental 
data of the electric and the magnetic properties can be explained in relation 
to the drastic change between the cluster and shell-model-like structures. 

In the following discussion, we consider two kinds of the fundamental 
effects of the cluster structures on the properties such as the electric and 
magnetic moments. One is caused by the spatial relative distance between 
clusters (spatial clustering effects) , and the other is concerned with the angular 
momentum coupling correlation of nucleons. As a typical example of the latter 
effect, we recall the so-called shell- model cluster in the SU 3 coupling shell- 
model [48,49] configuration. As we show below, it is found that the effect of 
the cluster structure on the magnetic dipole moments of Li and B isotopes is 
given only by the cluster coupling of angular momenta, while the quadrupole 
moments are effected by the spatial cluster as well as the cluster coupling of 
angular momenta. In order to extract the effect of the cluster coupling of 
angular momenta from the AMD wave functions, we have artificially made 
the inter-cluster relative distances in the AMD wave functions very small so as 
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to obtain the states in the shell-model limit. In the obtained shell-model limit 
states, the spatial cluster is not recognizable any more but only the effects of 
the cluster coupling of angular momenta persist. 



Table |VII| shows the results of [i and Q moments calculated with the spin- 
parity projected states from the shell-model limit wave functions, which are 
compared with the original AMD results. In the table we also list the ex- 
pectation values of the squared total-angular momenta J p , J n , J, the orbital- 
angular momenta L p , L n , L and the intrinsic spins S p , S n , S for protons, 
neutrons and for the total system, which have close relations with the spin 
configurations. We also calculate the z-components of the orbital-angular mo- 
menta and the intrinsic spins of protons and neutrons in the highest M states 
| JM) = |||). Since the /x moments in the shell- model limit are found to be 
almost the same as those of the original AMD, it is confirmed that the mag- 
netic dipole moments do not depend on the spatial clustering but are effected 
only by the cluster coupling of angular momenta. It is easily understood 
because the expectation values of linear terms of operators, J, L and S are 
mainly determined by the cluster coupling of angular momenta. As shown in 



the Table VII the magnitude of the total neutron intrinsic spin almost equals 
to because the intrinsic spins of the even neutrons are almost paired off. 
It means that the direct contribution to the [i moments from the neutrons 
is little. Then the \x moments of odd-even Li and B isotopes approximately 
consist of two terms from z-components of S p and L p , 

/i w n s + m (25) 

fi s = 5.58(2 ll'SWIm)/^' W = (2 2 1-^1 22)/^"' (26) 

where /in stands for the nuclear magneton. Figure p^ presents the components 
from the two terms jjl s and \xi in the total /i moments. 

In the Li isotopes, {i s is almost constant because |S P | = \ and S pz = 0.5 



(see Table VIII). The deviation of /i from the Schmidt value and its N 
dependence originate from the latter terms \i\ due to the orbital-angular mo- 
mentum of protons. The reason of the N dependence of \ii can be understood 
as follows. The total angular momenta |L| are always 1 and L z = 1.0 in Li iso- 
topes. However the state with the cluster coupling contains the components 
of |L n | 7^ states with non-zero total-orbital-angular momentum of neutrons, 
which makes L pz decrease in the highest |||) states. That is why 7 Li with 
the clustering structure has the smallest fi moments in odd-even Li isotopes. 
In other words, the N dependence of fi moments of Li isotopes is described 
by the non-zero L n of neutrons due to the clustering coupling effects. 

In the case of B isotopes, the deviations of the [i moments from the 
Schmidt value 3. 79// n for the p 3 / 2 proton orbit are not as small as the case of 
Li isotopes. In B isotopes the z-component (S pz ) has proved to be about 0.35 
which gives a smaller magnitude fj, s ~ 1.8/xat than the case of (S pz ) = 0.5. 
In the results in Table VII, it is found that the calculated values of (L 2 ) are 



almost constantly 2.8 in B isotopes and 2.0 in Li isotopes. It implies that 
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TABLE VII. The expectation values of the magnetic dipole operator fi, the electric quadrupole 
operator Q, the square of total orbital-angular-momenta of protons (L 2 ), the square of total or- 
bital-angular-momenta of neutrons (L 2 ), the square of total orbital-angular-momenta of the system 
(L 2 ), the total intrinsic spin of protons (S 2 ), the total intrinsic spin of neutrons (S 2 ), and the total 
intrinsic spin of the system (S 2 ) for the highest states |||) in the simplest AMD calculations of Li 
and B isotopes. Those in the shell-model limit are also listed. 
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2.0 


2.0 


0.00 


0.75 


0.75 




EXP. 


2.69 


40 














11 B 


AMD 


2.65 


34.0 


2.5 


3.6 


2.8 


0.04 


0.75 


0.78 




SML 


2.66 


25.9 


2.3 


3.4 


2.8 


0.04 


0.75 


0.78 




EXP. 


3.17 


37(4) 














13 B 


AMD 


3.17 


31.7 


0.0 


2.7 


2.7 


0.00 


0.75 


0.75 




SML 


3.18 


28.6 


0.0 


2.7 


2.7 


0.00 


0.75 


0.75 




EXP. 


2.66 


38(1) 














15 B 


AMD 


2.63 


34.3 


3.7 


3.8 


2.8 


0.00 


0.75 


0.75 




SML 


2.64 


22.5 


3.5 


3.7 


2.8 


0.00 


0.75 


0.75 




EXP. 


2.54 


— 














17 B 


AMD 


2.49 


42.2 


4.4 


4.1 


2.9 


0.07 


0.75 


0.81 




SML 


2.50 


22.6 


4.0 


3.7 


2.9 


0.33 


0.75 


0.77 




EXP. 


— 


— 














19 B 


AMD 


2.53 


50.8 


4.3 


4.2 


2.90 


0.00 


0.75 


0.75 




SML 


2.55 


24.9 


3.9 


3.8 


2.9 


0.00 


0.75 


0.75 


^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
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TABLE VIII. The expectation values of the z-components of the operators which effect on the 
jjl moments; the total orbital-angular-momenta of protons (L pz ), the total orbital-angular-momenta 
of neutrons {L nz ), the total intrinsic spin of protons (S pz ), and the total intrinsic spin of neutrons 

i m +■ In a cim t-v Incf A l\/l I I 

122/ 

Li and B isotopes. 



(S nz ), for the highest states |||) in the simplest AMD calculations and in the shell-model limit for 



^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 






\Spz) 


\&nz) 


(Lpz) 


\-L>nz) 


7 Li 


AMD 
SML 


0.50 
0.50 


0.00 
0.00 


0.36 

0.35 


0.64 
0.65 


9 Li 


AMD 
SML 


0.50 
0.50 


0.02 
0.02 


0.71 
0.72 


0.27 
0.26 


11 Li 


AMD 
SML 


0.50 
0.50 


0.00 
0.00 


1.00 
1.00 


0.00 
0.00 


11 B 


AMD 
SML 


0.34 
0.34 


0.00 
0.00 


0.74 
0.74 


0.41 
0.41 


13 B 


AMD 
SML 


0.37 
0.37 


0.00 
0.00 


1.13 
1.13 


0.00 
0.00 


15 B 


AMD 
SML 


0.33 
0.34 


0.00 
0.00 


0.77 
0.77 


0.40 
0.40 


17 B 


AMD 
SML 


0.33 
0.33 


0.00 
0.00 


0.68 
0.68 


0.49 

0.49 


19 B 


AMD 

SML 


0.33 
0.33 


0.00 
0.00 


0.73 
0.73 


0.45 
0.45 


^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
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FIG. 22. The effects due to the intrinsic spins [i s and of the orbital-angular momentum m on 
the total magnetic dipole moments [i of Li and B isotopes. The experimental data are also shown 
with square points. 
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FIG. 23. The Q moments of Li and B isotopes for the original AMD wave functions and the 
shell-model limit described in the text. Squares are the experimental data. The spatial clustering 
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the original results (open circles). 
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the wave functions of Li isotopes are the almost pure L = 1 states, while 
those of B isotopes contains L = 1 states and L = 2 whose components are 
easily estimated to be about 20%. The mixing of the L = 2 states in B iso- 
topes makes S p tilted from the z direction in the \JM) = |||) states as to 
S pz ~ 0.35. The reason of the pure L = 1 state in Li and the mixing ratio of 



the L = 2 states in B are described in Ref. |3C] in more detail. The second 
term [i[ causes the N dependence of [i moments of B isotopes as similar way 
with Li isotopes. Even though the magnitude |L| and the z-component (L z ) 
of the total-orbital-angular momentum is approximately constant in the B 
isotopes, the clustering gives the components of L n ^ states which reduce 
the magnitude of {L pz ). Because of such a effect of the clustering coupling 
on the \i\ term, the magnetic moments of 11 B, 15 B, 17 B and 19 B are smaller 
than the one of 13 B which has the shell-model like structure. Since the closed 
neutron shell in the nuclei 13 B with N = 8 make S n = L n = J n = 0, the 
\x moments of 13 B is the largest in B isotopes. Similarly to Li isotopes, it is 
concluded that the N dependence of fj, moments of B isotopes is understood 
by the cluster coupling effect. 

In contrast to the magnetic dipole moments, the electric quadrupole mo- 
ments are sensitive to the relative distance between clusters. For the sake of 
estimating the spatial cluster effects on Q moments, we try to decompose the 
calculated Q moments into two components: The first component is originated 
by the spatial clustering, and the other is due to other properties including 
the cluster coupling of angular momenta. We consider that the second com- 
ponents are given by the Q moments calculated by the shell-model-limit wave 
function where the spatial clustering has been already removed. They are 



shown in Table VII together with the Q-moments of the AMD calculations. 
We display the two components in Q moments in the Fig. p3, which is helpful 
to analyze the N dependence of the each component. As for the B isotopes 
the calculated Q moments in the shell model limit are 25.9, 28.6, 22.5, 22.6 
and 24.9 e-mb for 11 B, 13 B, 15 B, 17 B and 19 B, respectively. It shows that 
the second components in these nuclei except for 13 B are smaller than the 
one in 13 B with A = 8. With similar argument as for /^-moments, we expect 
that the reduction of the Q-moments in other B isotopes than 13 B should be 
explained by the mixing of the component of the non-zero neutron orbital- 
angular momentum. By subtracting these second components from the total 
Q-moments ( namely the Q-moments of the AMD calculations ) , we can esti- 
mate the contribution of the first component due to the spatial clustering as 
8.1, 3.1, 11.8, 19.6 and 25.9 e-mb for n B, 13 B, 15 B, 17 B and 19 B, respectively. 
This component is smallest in 13 B and becomes larger as the neutron number 
increases toward the neutron drip-line. Such dependence on A of the first 
component is indeed consistent with the clustering development mentioned 
in the previous subsection. In the case of Li isotopes the first component is 
largest in 7 Li and decreases toward n Li as the clustering structure weakens. 
Thus it is proved that the systematic A-dependent features of experimen- 
tal data for the electric magnetic properties are quantitatively explained by 
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the structure change given by our AMD results. The reader is referred to Ref. 



1 30|| for the detailed discussions. 



3. Opposite deformation between protons and neutrons 

As mentioned above, AMD results suggest that the opposite deformations 
between protons and neutrons may be found in the proton-rich C isotopes. In 
Fig. ^4|, we illustrate the deformation parameters of protons and neutrons in 
proton-rich C isotopes. It is notable that the neutrons prefer the prolate and 
the triaxial deformations rather than the oblate shape in these nuclei because 
of the neutron number N = 3,4,5. On the other hands, the protons prefer 
oblate shape in C isotopes. As a result the disagreement between the proton 
and the neutron deformations is found in proton-rich C. Our purpose here is 
to confirm the disagreement between the proton and the neutron shapes by 
the help of the electric quadrupole moments and transitions in C and the ones 
in the mirror nuclei. 

First we discuss Q moments of n C and the mirror nucleus n B. Based on 
mirror symmetry for proton and neutron deformations, we compare the proton 
and the neutron deformations by analyzing the ratio of the electric quadrupole 
moment Q in n C to the one in 11 B. We introduce the well-known approximate 
relation between the electric quadrupole moment Q in the laboratory frame 
and the intrinsic quadrupole moment Qq; 

3K*-J(J + 1) 
Q ~ g °(2J + 3)(J + l)- (27) 

By using Eq. Bl] we can express the intrinsic electric quadrupole moment Qq 
as follows in the first order of the deformation parameter f3 p , 



Qo = \I^^Ze(3 p cos lp Ri (28) 

where f3 p and "f p are the deformation parameters for the proton density, and 
Z and R e are the proton number and the charge radius, respectively. In- 
stead of the usual deformation parameter /3 p in the usual equation Qo = 

J ^--^ZeftpRz, we introduce the effective deformation parameter f3 p cos^ p in 
Eqj2^ which is necessary for the system with different proton and neutron 
shapes as described below. We explain the appropriate principal axes in the 
nucleus where the different proton and neutron shapes coexist as shown in 
schematic Figure g{|. For example in the nucleus with oblate proton and 
prolate neutron deformations, the approximate symmetry axis x for protons 
usually differs from the approximate symmetry axis z for neutrons so as to 
make the largest overlap between the proton density and the neutron den- 
sity. In many cases, it has been found that a symmetry axis z for the prolate 
neutron density should be chosen as the principal axis Z of the total intrinsic 
system for the total-angular momentum projection in generating the lowest 
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FIG. 24. Deformation parameters (/3, 7) for proton and neutron densities in the intrinsic states 
of proton-rich C isotopes with A = 9 ~ 11 are given by triangles and squares, respectively. The 
mass number A is indicated beside each point in the figure. The interaction (b) is used. 




J^ state with an approximately good K quantum. In such cases, the usual 
formula for Qq should be modified by using the effective deformation to z-axis, 
/3 p cos7p instead of [3 P . In other words, an oblate deformation gives a smaller 
contribution to the intrinsic quadrupole moment Qo than is expected. 

Assuming these simple approximations the ratio of the Q moment in n C 
to that in n B is represented by the product of three terms, the ratio of proton 
numbers, that of the proton deformation parameters, and the one of the charge 
radii. When we assume the mirror symmetry for the deformation parameters 
and replace the deformation parameter /3 p cos7 p ( 11 B) of the proton density in 
11 B by /3 n cos7 n ( 11 C) of the neutron density in 11 C, the ratio of Q is written 
as, 



QC l c) 

Q( n B) 



Z = 6 
Z = 5 



/ /? p cos7 P ( n C) \ ( R e ( ll C) 

U n COS7n( n C)i li?e( n B) 



We take the third term for the charge radii to be unity since n C and n B 
are the nuclei close to the stability line. If the neutron deformation agrees 
to the proton deformation in 11 C, the second term gives no contribution to 
the ratio of Q moments and the ratio can be explained only by the charge 
ratio 1.2. However the ratio 0.88 deduced from the experimental data is less 
than a unity and is inconsistent with the Z ratio because the experimental 
data of Q( 11 C)=34.8 e-mb is smaller than Q( 11 B)=40.7 e-mb. According to 



(29) 
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FIG. 25. Schematic figures of the nucleus with oblate proton and prolate neutron defor- 
mations. In the body-fixed flame x, y, z is chosen so that moments of inertia obey the relation 

T < T < T 





oblate 



prolate 



AMD calculations, this problem can be resolved by taking into account the 
difference between the intrinsic deformations of proton and neutron densities 
in n C. Through the second term in Eq.(f2S|), the ratio of Q moments reflects 
the difference of intrinsic shapes between proton and neutron densities. As 



shown above in Fig. 24 the proton deformation is oblate while the neutron 
density becomes a triaxial shape in 11 C. Since /3pCOS7p( 1:L C) is smaller than 
/3 n cos7 n ( 11 C), the second term in Eq. 29 becomes less than unity, which 



cancels the effect of the first term of the Z ratio. 

We make further quantitative discussion by the use of theoretical values 
of deformation parameters in the intrinsic states obtained with AMD. In the 
present results, the ground state of n C with J = 3/2 is obtained by a total- 
angular-momentum projection on a state \J = 3/2, K = 3/2} with respects 
to the principal z axis with the minimum moment of inertia. Using the 
theoretical values of f3 p cosj p and /? n cos7 n shown in Fig.^ we can estimate 
the ratio of Q moments with the first and the second terms in Eq.[29J. The 
estimated ratio is found to be 0.87 which is as small as the value of 0.88 
deduced from the experimental data. In fact the theoretical results of electric 
Q moments for total-angular momentum projected states are 20 e-mb for 
n C and 34 e-mb for 11 B, which are consistent with the experimental data of 
Q( n C)< Q( n B) (Table.©. 

It is concluded that the difference of the intrinsic deformation of the proton 
density in n C from that in n B (an oblate shape in n C and a triaxial shape in 
n B) makes significant effects to the ratio of the Q moment of n C to that of 
11 B. By assuming the mirror symmetry, it is theoretically suggested that the 
disagreement between proton and neutron deformations in n C is supported 
by the experimental fact of the ratio Q( n C)/Q( n B)/ less than a unity. 

Next we make similar analysis of the deformations for 10 C. We find that 
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TABLE IX. Electric quadrupole moments and transitions of proton-rich C isotopes and the 
mirror nuclei. Calculations are with MV1 force (m = 0.576) and the experimental data are taken 
fromRef. I6|. 



nucleus 




level 


electric 
exp. 


Q 


moments 

theory 


li c 
n B 




3/2" 
3/2~ 


34.3 e-mb 

40.7(3) e-mb 




20 e-mb 
34 e-mb 


io c 
10 Be 




2+ 
2+ 


I 




—38 e-mb 
—65 e-mb 


9 C 
9 Li 




3/2- 
3/2" 


-27.8 e-mb 




-28 e-mb 
—27 e-mb 


nucleus 




level 


E2 transition strength 
exp. theory 


li c 
n B 


5/2 
5/2 


" -> 3/2- 
" -> 3/2- 


13.9(3.4) e 2 fm 4 




6.8 e 2 fm 4 
11.3 e 2 fm 4 


io c 
10 Be 


2+^0+ 
2+^0+ 


12.3(2.0) e 2 fm 4 
10.5(1.0) e 2 fm 4 




5.3 e 2 fm 4 
9.5 e 2 fm 4 


9 C 

9 Li 


1/2 
1/2 


" -» 3/2" 
" -> 3/2- 


- 




5.7 e 2 fm 4 
7.2 e 2 fm 4 



the difference between proton and neutron deformations in U C is important 
to understand the ratio of the El transition strength B(E2; 2 + —>■ + ) in 10 C 
to that in the mirror nucleus 10 Be. Assuming mirror symmetry, the ratio of 
B(E2) is approximated similarly to Eq.^ as, 

2 



B(E2; W C) /Z = 6 
B(E2; 10 Be) ~ \Z = 4 



/ /? p cos 7p ( 10 C) 
U n cos 7 „( 10 C) 



Re( w C) 
R e ( l °Be) 



The first term of the charge ratio (6/4) 2 =2.25 is much larger than the ratio 
of experimental values 12.3(2.0)e 2 fm 4 /10.5(1.0)e 2 fm 4 =1.2(0.3). The reason 
why the square of the charge ratio fails to reproduce the ratio of B(E2) in 
the mirror nuclei 10 C and 10 Be is because of the disagreement between proton 
and neutron deformations in 10 C. 

In the intrinsic state of 10 C, the proton density deforms oblately with 
PpCosjp = 0.28 while the neutron deformation is prolate with a larger value of 
the effective deformation parameter /3 n cos7 n =0.49 (FigJ24|), which makes the 



second term in Eq.3C less than unity. If the third term is omitted, the ratio 
is roughly estimated as, 

B( ; E \ C) , = 2.25 x (0.28/0.49) 2 ~ 0.75. 
B(E2; l0 Be) v ' ' 

Although the ratio estimated above is smaller than the ratio 1.2 deduced from 
the experimental data, it is found that the reduction of the ratio of the proton 



(30) 



(31) 
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numbers is made by the ratio of the deformation parameters The microscopic 



calculations of Q moments with AMD are shown in Table |IX| and are compared 
with the experimental data. Since the calculations underestimate the value of 
B(E2; 10 C), the ratio B (E2; 10 C) / B (E2; w Be) is smaller than a unity. In the 
results with VAP calculation which will be mentioned later, the theoretical 
value of B(E2; 10 C) is improved. 

As for the 9 C and the mirror nucleus 9 Li, the present prediction is that 
the Q moments of 9 C is smaller to the one of 9 Li. We should note that the 
present results do not include the effect of the long tail of valence nucleons. 
However, if the proton halo of 9 C exists because the separation energy of 
protons is small, the orbits of the valence protons may give the large effect on 
the Q moments. 

4- Neutron skin and halo 
The presence of a neutron skin structure has been discussed for a long 



time. Recently thick neutron skins have been reported in He isotopes [50] 
and in 20 N [ET] by the help of the experimental data of interaction cross 
sections. The appearance of the neutron skin has been also shown in the 
comparison of neutron radii with proton radii along a chain of Na isotopes 
by combining the data of the isotope-shift for charge radii with those of the 
matter radii deduced from interaction cross sections |p2j] . Also the theoretical 



studies of skin structures have been tried in unstable nuclei p2}3C,44,53|. 

The present results suggest that in the neutron-rich nuclei of B and C 
isotopes the density of neutrons stretches far toward the outer region. The 
simple AMD calculations predict the presence of "neutron skin structure", 
which is the surface region with rather high neutron density but low proton 
density. In particular, C isotopes are expected to have thicker skins than 
those of B isotopes because the neutron-rich C have no developed clustering 
structure as is seen in the proton density of 20 C which remains in the inner 
region as compact as that of stable C nuclei. One of characteristics of C iso- 
topes is the stationary structure of protons in spite of the drastic change of 
neutron structure along the increase of neutron number. It is in contrast with 
neutron-rich B isotopes which are predicted to have the clustering structure. 
According to simple AMD calculations, in a series of A = 14 isotones the 
neutron skin which is developed in 20 C weakens with the increase of the pro- 
ton number because the mean-field for the neutrons given by proton density 
becomes deeper to decrease the neutron radii. 

Figure |2^ shows the proton and the neutron densities of C isotopes as a 
function of radius. The difference between proton and neutron densities in 
the surface region around r ~ 4 fm enhances as the neutron number increases 
from N = 10 toward the neutron-drip line. In the surface region of 20 C, the 
line for the neutron density (solid) seems to be shifted outward by about 1 
fm compared with the line for proton density (dashed). 
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FIG. 26. The densities p p of protons and p n of neutrons as a function of radius in C isotopes. 
The densities of the + states are obtained with AMD calculations by using MV1 force with 
m = 0.63. The proton density of 12 C is written by the dotted lines. The densities of in 19 B are 
also shown together with the proton density of 13 B. 
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FIG. 27. Theoretical values of root-mean-square radii of protons (dotted lines) and neutrons 

(solid lines). Calculations are the results with the interaction (b) MV1 (m = 0.576) and with the 

MV1 force with m = 0.63. 
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As already mentioned the recently measured radii of C isotopes [43] agree 
well systematically with the simple AMD calculations except for 19 C which 



is expected to have a neutron halo structure (Fig. |10|). We show the root- 
mean-square radii of neutrons and protons separately in Fig.^7]. It is found 
that the neutron radii become larger and larger in the region heavier than 
14 C while the proton radii are rather stable with the increase of the neutron 
number. Although the radii depend on the adopted interaction parameters, 
in both calculations with Majorana parameter m = 0.576 and m = 0.63 the 
difference between the proton radius and the neutron radius in 20 C is more 
than 0.3 fm. According to the present results, the increase of the matter 
radii in the neutron-rich C isotopes is mainly due to the development of the 
neutron skin structure. 

The radial behavior of the densities of protons and neutrons is related 
closely with the single-particle energies. We calculate the single-particle ener- 
gies and the single-particle wave functions by diagonalizing the single-particle 
Hamiltonian with the analogy to Hartree-Fock theory. First we transform 
the set of single-particle wave functions </?j of the solved of AMD wave func- 
tion into an orthonormal base ip a . The single-particle Hamiltonian can be 



constructed by the use of the orthonormal base as follows [44,37,3 



7 

1 A 
+ r ^2( l Pa l P~/ l P8\v3\ l Pf3 l P~f l P5 + fig'Pp'P'y + Vi&S&P ~ VpV&V-i ~ l P~f l Pf3 l P6 ~ 'PS&'y'Pp), (32) 

7,(5 
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where the Hamiltonian operator is written by a sum of the kinetic term, 
the two-body interaction term and the three-body interaction term; H = 
J2i t + J2i<j $2 + J2i<j<k ^3- We note that the single-particle energies defined 
as above can be positive because the model space for single particle wave 
functions is restricted to the Gaussian wave packet giving rise to an artificial 
wall due to the zero-point kinetic energy of the packet which prevents the 
single-particle wave function from escaping out of the nucleus. In Fig.p3 we 
present single-particle energies in C isotopes which are calculated from the 
obtained intrinsic AMD wave functions. It is found that in the neutron-rich 
C isotopes heavier than 14 C the Fermi energy is at most a few MeV. In the 
nuclei near the neutron-drip line, many valence neutrons occupy the higher 
orbits a few MeV below the zero energy which correspond to sd orbits. These 
weakly bound neutrons build the neutron skin structure. In the even-odd 
neutron-rich C isotopes such as 19 C, the energy of the last valence neutron 
is about zero energy. The possible existence of halo structure in 19 C is sug- 
gested experimentally by the measurements of the longitudinal-momentum 
distribution of 18 C after the one-neutron breakup of 19 C [13] and also by 



the interaction cross sections [43|. Although the halo structure may not be 
seen in the present results because the AMD wavefunction is not sufficient to 
describe a long tail of the halo tail, considering the small binding energy of 
valence neutrons we naturally expect that the halo structure may appear if 
the single-particle wave function of such a loosely bound valence neutron is 
described more precisely than the present model space. 

On the other hands, the protons are deeply bound in the neutron-rich 
C. The binding energies of protons grow rapidly from 9 C to 14 C because the 
number of neutrons in p shell increases. In the heavier region from 14 C to the 
neutron-drip line, the proton energy becomes deeper and deeper gradually. 
Taking into account the potential depth of protons and the enlargement of 
the neutron density distribution, the kink of radii at 14 C seems to be a natural 
phenomenon which reflects the shell effect of neutrons at N = 8. 

V. STUDY WITH VAP 

Recently the structures of the excited states as well as the ground states 
are very attractive in the study of unstable nuclei. It is natural that various 
molecule-like states may appear in the excited states of light unstable nuclei 
because the excitation due to the relative motion between clusters is important 
in the light nuclear region. We apply the method of the variational after spin- 
parity projection in the framework of AMD to the light unstable nuclei for 
the aim to make systematic study of the structure change with the increase 
of the excitation energy. The formulation has been already explained in the 
section |fl]. The applicability of the framework for the study of the excited 
states in the stable nuclei has been confirmed in Ref. |38| on the structure of 
12 C. Here we study the structures of the excited states of 10 Be and 12 Be. 
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FIG. 28. The single-particle energies of the intrinsic states for the normal parity states of C 
isotopes. The method to calculate the single-particle energies is explained in the text. The energies 
of proton (neutron) orbits are presented on the left(right). The MV1 force with m = 0.63 is used. 
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A. 10 Be 

10 Be, one of the challenges in the study of light unstable nuclei, has been 
investigated experimentally not only by use of the unstable nuclear beams but 
also in such experiments as the transfer and pick-up reactions. The recent 
experiments of the charge exchange reactions 10 B( 3 He,t) 10 Be [J54J let us know 
the strength of the Gamow- Teller transitions to the excited states of 10 Be. 
These new data of (3 transition strength which are deduced from the cross 
sections at the forward angle are very helpful to study the structure of the 
excited states. The structures of 10 Be have been studied hard also theoret- 
ically by microscopic calculations, for example, shell models |55,p6|, cluster 



models [OJlE.47], Hartree-Fock [|57|], and antisymmetrized molecular dynam- 
2S,30,37|. We study the structure of the excited states of 10 Be with the 



1CS 



VAP calculations in the framework of AMD. 



1. Results 



In this subsection we display the theoretical results of the excitation en- 
ergies, E2, El, and /3 transitions which can be directly compared with the 
experimental data. The detail of the structures is discussed in the next sub- 
section. 

The adopted interaction for the central force is the case 3 of MV1 force 



[41]. The adopted parameters are m = 0.62, b = h = for the Majorana, 
Bertlett and Heisenberg terms of the central force and the strength of the 
spin-orbit force ui = — un = 3000 MeV (interaction (g)). Trying another set 
(h) of parameters with m = 0.65, b = h = and ui = —uu = 3700 MeV, 
we did not find significant differences in the results. The set of parameters of 
case(g) is the one adopted in the work on 12 C []38|| . On the other hands, the 
VAP calculations with the set of interactions case(h) reproduce the abnormal 
parity of the ground state of 11 Be. The optimum width parameters v of wave 
packets are chosen to be 0.17 fm -2 for case(g) and 0.19 fm ~ 2 for case(h) 
which give the minimum energies in VBP calculations of 10 Be. The resonance 
states are treated within a bound state approximation by situating an artificial 
barrier out of the surface. 

The lowest J^ states are obtained by VAP calculations for P mk ,$amd 
with {J±,K')= (0 + ,0), (2+,0), (3+,+2), (4+,0), (1",-1), (2",-l), (3~,-l), 
(4~, —1). Considering O-j" state to be + state in the second K w = 0^ band, 
the O2" state is calculated by VAP as the higher excited state orthogonal to 
the lowest Of state as explained in the subsection |IIF , That is to say that 



the 0^~ state is obtained by VAP for $(Z) in Eq. || with (J^if',™) 
(0 + ,0, 2). In the case of higher 2 + states, we impose the constraint that the 
approximately principal z-axis of the intrinsic deformation equals to the 3-axis 
of the Euler angle in the total-spin projection. According to VBP calculations 
the second 2 + state is described as the band head of the lowest K w = 2 + band. 
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Therefore we construct the state 2^~ by choosing (J ,K') of P MK i^amd as 
(J^,K') = (2 + ,+2) under the constraint on the principal z-axis. Because 
of the constraint and the choosing of the K' quantum number it keeps the 
approximately orthogonality to the lowest 2 + state with (J ,K') = (2 + ,0). 
The third 2^ state is easily conjectured to be 2 + state in the second K w = Of 



band like O2 state. We obtain the 2J state by VAP for 3>(Z) in Eq.[L5] with 
(J^, K', n) = (2 + , 0, 2) by imposing the constraint so as to make the principal 
z-axis equal with the 3-axis in spin projection. It means that the orthogonal 
condition of 2g~ to 2f is kept by superposing two wave functions as described 



in the subsection II F, while the orthogonality to 2^ (K 77 = 2 + ) is taken into 
account by choosing the different K quantum K' = 0. 

The binding energy obtained with the case(g) interactions is 61.1 MeV, 
and the one with case(h) is 61.3 MeV. The excitation energies of the results 
are displayed in FigJ2^. By diagonalization of the Hamiltonian matrix the 
excited states 4^, 6^ are found in the rotational band K n = 0^~, and 5~ state 
is seen in the K 11 = 1~ band. Comparing with the experimental data, the 
level structure is well reproduced by theory. Although it is difficult to estimate 
the width of resonance within the present framework, the theoretical results 
suggest the existence of 3 + , 4 + , 6 + and 5~ states which are not experimentally 
identified yet. The excited levels can be roughly classified as the rotational 
bands K* = Of, 2+, 0^ and 1~ which consist of (Of, 2+, 4+), (2+, 3^), (0^, 
2g, &21 ®i) an< i (I - ) 2 _ , 3 _ , 4 _ , 5 _ ), respectively. The intrinsic structures 
of these rotational bands are discussed in detail in the next section. 

The data of transition strength are of great help to investigate the struc- 
tures of the excited states. The results with the interaction case(g) and the 
experimental data of El and El transition strength are listed in Table [X|. 
The theoretical values well agree with the experimental data. The strength 
B{E2) for 10 C;2^ — ► 0^ is simply calculated by the wave function of 10 C sup- 
posed to be mirror symmetric with 10 Be. The present result for E2 strength 
of 10 C;2^ — ► 0^ is better than the work with simple AMD calculations (see 



Fig. IX) [45]. As for the values with a shell model, (0 + 2)Tiuj shell model cal- 



culations from the reference |5£] with effective charges e n = 1.05e, e u = 0.05e 



are also listed. The shell model calculations qualitatively reproduce some 
experimental data of the E2 properties of low-lying levels. 

The strength of the (3 decays of Gamow-Teller(GT) type transitions can 
be deduced from the cross sections at the 0° forward angle of the charge 
exchange reactions which have been measured recently [E3]. These new data 
for the Gamow- Teller type j3 transitions are very useful probes to discuss the 



structures of the excited states of unstable nuclei. Table |X| shows the values of 
B{GT). The experimental values for the j3 transitions from 10 B(3+) to 10 Be* 
are deduced from the data of the reaction 10 B(t, 3 He) 10 Be. In the present 
calculations, the wave functions for the neighbor nucleus 10 B are calculated 
with VAP where (J ,K') is chosen to be (3 + , —3) for the ground 3^ state and 
(1 + , —1) for the if state. 10 Be and 10 B are calculated with the case (g) and 
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FIG. 29. Excitation energies of the levels in 10 Be. Theoretical results are calculated by the 
diagonalization of the states obtained with VAP by using the interaction case (g) . 
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TABLE X. i?2 and El transition strength. The theoretical results of VAP calculations with 



the interaction case (g) are compared with the experimental data [58]. The shell model calculations 
are quoted from the work with the (0 + 2)hu> shell model in the reference [|5S[|. 



transitions 


Mult. 


exp. 


present VAP 


shell model 


10 Be;2+ -> 0+ 
10 Be;0+ -* 2+ 
10 Be;0^ -► li 


£2 
£2 
El 


10.5±1.1 (e fm 2 ) 

3.3±2.0 (e fm 2 ) 

1.3±0.6xl0" 2 (efm) 


11 (e fm 2 ) 

0.6 (e fm 2 ) 

0.6xl0" 2 (efm) 


16.26 (e fm 2 ) 
7.20 (e fm 2 ) 


10 C;2+ -> 0+ 


£2 


12.3±2.0 (e fm 2 ) 


9 (e fm 2 ) 


15.22 (e fm 2 ) 
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TABLE XL B{GT) values of (3 decays which are the square of the expectation values of 
Gamow- Teller operator. The experimental data are the values'"' deduced from the cross sections 
of 10 B(t, 3 He) 10 Be* at 0° forward angle |J4J and the one' ' from |60| . The theoretical results are 
obtained with case (g) and case (h) interactions for 10 Be and 10 B. 





exp. 




initial 


final 




{J\E X ) (MeV) 


{J\E X ) (MeV) 


B(GT) 


10 B(3+,0) 


10 Be(2+,3.37) 


0.08± 0.03 a ) 


10 B(3+,0) 


10 Be(2^,5.96) 


0.95± 0.13 a ) 


10 B(3+,0) 


10 Be(2+ or 3+, 9.4) 


0.31± 0.08") 


10 C(0+,0) 


10 B(1+, 0.72) 


3.44 fe ) 


theory case(g) 


10 B(3+) 


10 Be(2+) 


0.02 




10 Be(2+) 


1.1 




10 Be(3+) 


0.40 




10 Be(4+) 


0.08 




10 Be(2 3 +) 


0.03 


10 Be(0^) 


10 B(1+) 


2.9 


theory case(h) 


10 B(3+) 


10 Be(2+) 


0.00 




10 Be(2+) 


0.92 




10 Be(3^) 


0.38 




10 Be(4+) 


0.10 




10 Be(2 3 + ) 


0.00 


10 Be(0+) 


10 B(1+) 


2.5 



(h) interactions. The theoretical values reasonably match to the experimental 
data. Since the data for 10 B(3 + ) — > 10 Be(9.4MeV) well correspond to the 
theoretical value of 10 B(3+) -> 10 Be(3i"), the excited level of 10 Be at 9.4MeV 
is considered to be the 3f state. The strength of these GT transitions from 
10 B(3 + ) is governed by the configuration of the ground state of 10 B which is 
understood as the state 3 + with \K\ = 3 in p-shell in the simple shell model 
limit. It is natural that the transitions to 2^ and 3 J" states in the K w = 2 + 
bands of 10 Be are strong while the transitions to the states in the K n = + 
bands are weaker. That is why the predicted B(GT) for 10 B(3 + ) ^ 10 Be(23~) 
is small because 10 Be(23~) is the state in the K w = 0^ band constructed by 
the linear structure with the developed clustering. The result of B(GT) for 
10 Be(0^~) — > 10 B(1 + ) is consistent with the experimental value of the (3 decay 
from the mirror nucleus 10 C(0^) — > 10 B(1 + ). 
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2. Intrinsic structure 

Here we discuss the structure of the excited states by analyzing the wave 
functions. Even though the state calculated by VAP mixes with each other af- 
ter the diagonalization of the Hamiltonian matrix, the state P MK ,^AMD{'^ri iz ) 
projected from a Slater determinant obtained in VAP with (J ,K',n) is the 
main component in the final result of the J„ state. In this section we consider 
the Slater determinant $amd(Z^ ± ) as the intrinsic state for the J„ state. 

In the excited states, various kinds of structures are found. Here we 
analyze the structures of the intrinsic states $imd(Z„ )• It is found that 
the excited levels are classified into rotational bands as 0J~, 2f, 4J~ states in 
K n = Of band, 2^, 3f in K* = 2+ band, 0^, 2^, if, 6f in the second 
K w = Of band and 1", 2", 3", 4~, 5" in K n = l~ band. Particularly the 
molecule-like states with the well-developed 2a cores construct the rotational 
bands i^ = OJ and 1~ in which the level spacing is small because of the large 
moments of inertia. The density distributions of matter, protons and neutrons 
in the intrinsic states $AMfl(Z^) are presented in Fig. |3(| We found the 
2a + 2n structures in most of the intrinsic states. The density of protons 
indicates that the clustering structure develops more largely in 1~ than 0^~ 



and most remarkably in Of ■ As seen in Fig. 30, the intrinsic structure of Of 
state has an axial symmetric linear shape with the largest deformation, while 
1~ state has an axial asymmetric shape because of the valence neutrons. The 
structures of 0^ and V[ states are similar to the ones of the previous work 
with the simplest version of AMD |4J]. The increase of the degree of the 
deformation along 0^, 1~ and 0^ is consistent with the previous works such 
as [18,^7,3^]. In the K = 1~ band, the deformation toward the prolate 
shape shrinks as the total spin J increases. 

In the K 71 = 0^ band, the development of the 2a cores weaken with 
the increase of the total spin due to the spin-orbit force. The reduction of 
the clustering structure is more rapid in the case of interaction (h) with the 
stronger spin-orbit force. As a result, the 2^" and 4^" states in the case (h) 
interactions contain the dissociation of a. Regarding the dissociation of the 
a cores, the structures of those states are sensitive to the strength of the 
spin-orbit force. 

B. Behavior of valence neutrons 

Even though we did not assume the existence of any clusters in the model, 
we have found that the most of the intrinsic states of 10 Be contain 2a + 2n 
structures. We study the behavior of the valence neutrons surrounding 2a 
by analyzing the single-particle wave functions to understand the role of the 
valence neutrons in the neutron-rich Be nuclei. Considering that an intrinsic 
state is approximately written by a Slater determinant &amd, the single- 
particle wave functions and the single-particle energies of an intrinsic state 



08 



p/2 p p 



•n 




•• • 



• • 


•• 


f 





? + 



21 



Oo + 



1 



3" 



10 fm 10 fm 10 fm 

FIG. 30. The intrinsic structure of the excited states of 10 Be obtained by VAP calculations. 
The density distribution of matter, protons and neutrons of the intrinsic states are shown at left, 
middle and right, respectively. The density is integrated along the axis perpendicular to adequate 
planes. The figures are for the results with the interaction (g). 



• • 




• • 




•% 




59 



are determined by diagonalizing the single-particle Hamiltonian which has 



been already explained in IV B 4. 



The single-particle energies in the 0^ , 1~ and Og states are shown in 
Fig. [n|. In each state four neutron orbits from the bottom correspond to 
those for the neutrons in the 2a clusters. The level spacing of these four 
lower orbits becomes smaller in l - than in 0^ and smallest in the 0^ state 
with the increase of the distance between clusters. We consider the last two 
neutrons in the higher orbits as the valence neutrons surrounding 2a cores. 
We display the density distributions of the single-particle wave functions for 



the two valence neutrons in the left column of Fig. 32. Figures in the middle 



and right columns of Fig. |32j are for the normalized density of the positive 
and the negative parity components projected from the single-particle wave 
functions, respectively. By analyzing the single-particle wave functions it is 
found that two valence neutrons of the 0^~ states contain the negative parity 
components more than 80% (the right column in figure ^) which seem to be 
it bonds(Fig. ^(a)) in terms of molecular orbits. On the other hands in the 
case of Oj" band, the last 2 neutrons are predominantly in the positive parity 



orbits, which are analogous to the a bonds (Fig. 33 (b)). In the 1~ band, 
each valence neutron contains both the positive parity component like a and 
the negative parity one similar to n. Since the parity of the total system is 
negative in the 1~ band, the states are considered to have one neutron in a 
orbit and the other neutron in n orbit after the parity projection. Therefore 
roughly speaking, the 0^, 1~ and 0^~ states are understood as 2a and two 
valence neutrons in it 2 , air and a 1 orbits, respectively. The interesting point 
is that the valence neutrons play important roles to develop the clustering 
structure in the excited bands. We can argue that the clustering develops 
in the 1~ band and mostly in 0^~ band owing to the a orbits of the valence 
neutrons, because the a orbit prefers the prolately deformed system as to gain 
its kinetic energy. This idea originates from the application of the two centers 



shell model to the 2a dimer model of Be isotopes by W. von Oertzen |61] and 
consistent with the argument in the work with the method of AMD+HF by 
Dote et al. §7|. 

From the view points of the shell model, the levels J n =0^, 2^~, 4^ in the 
K n = Of band and 2^, 3^" in the K n = 2f are the OHlo states. These states 
are dominated by the components with the same [442] spatial symmetry in 
the supermultiplet limit. On the other hands, the negative parity states in the 
K w = 1~ are dominated by the lhto configurations as for the neutron p-shell. 
In the excited states in the K n = 0^" band, the 2hu configurations with 2 
neutrons excited are prime. It is consistent with the rotational band seen in 



the 2%oj states of shell model calculations [55]. However the well-developed 
clustering structures in the K 77 = 1~ and the K w = 0^ bands contain con- 
siderably the higher configurations which are not included explicitly in the 
wavefunction of the Vhuj model space and (0 + 2)Thjo model space of the shell 



model. Also in the Skyrme Hartree-Fock calculations 57], the Oj state is 
found to be a state with 2 neutrons in sd-orbits. 
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FIG. 31. Single-particle energies in the intrinsic systems of the 0J , 1 and 0<J~ states. The 
energies of protons(neutrons) are displayed in the left(right) side in each figure. 

C. Results of 12 Be 



In 12 Be, one of the interesting problem is the vanishing of the magic 
number N = 8 which is closely related with the abnormal parity of the ground 
state of the 11 Be. Another subject on 12 Be is the molecule-like structure in 



the excited states which is reported in Ref. [28|. 

In present calculations, we adopt the interaction (h) with which we can 
reproduce the spin parity of the ground state l/2 + of n Be. For comparison, 
we also use the interaction (g) which contains the weaker spin-orbit force than 
the interaction (h) and fails to reproduce the parity inversion in 11 Be. In the 
theoretical results of the energy levels we found many rotational bands such 

O2". We present the density distributions 



as K m 



07, K* 



>i ■ 



17 and K* 



of the intrinsic structures in Fig. 34 for the results with the interaction (g). 
It is surprising that the intrinsic state of the ground band K 71 = Of is the 
largely deformed molecule-like state with the dominant 2p-2h state in terms 
of neutron p-shell which is analogous to the second + state in 10 Be, while 
the second band K 71 = 0^ is dominated by the normal Ohco state with the 
rather spherical shape due to the neutron closed shell. It means the abnormal 
phenomena that the s<i-orbits energetically intrude in the p-orbits and the 
magic number of N = 8 vanishes in the 12 Be. On the other hands, adopting 
the interaction (g) we obtain the rather spherical ground state with the closed 
neutron p-shell, and the deformed structure with 2p-2h appears in the second 
+ state lying a few MeV above the ground state. Such a difference in results 
between the adopted interactions (g) and (h) originates from the strength of 
spin-orbit force. The stronger spin-orbit force in (h) prefers the closed P3/2- 
shell and breaks the p-shell in Be isotopes. With both interaction (g) and 
(h) we find the K n = 1~ band with the developed clustering structure in the 
low-lying levels. The K w = + band with the 2p-2h configurations reaches 
to the J 77 = 8 + state, and the band terminal of K w = 1~ is found to be a 
J 71 = 6~ state. The spins 8 + and 6~ at the band terminal are the same as 
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FIG. 32. Density distribution of single-particle wave functions of the valence two neutrons 



in the intrinsic states of 0^ , 1 



and Oj (left column). The method to extract the single-particle 



wave functions are explained in the text. The middle and right columns are for the density of the 
positive parity and negative parity components projected from the single-particle wave functions, 
respectively. The wave functions projected into the parity eigen states are normalized in the 
presentation. 
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(b) 



(5^*C3^«c3 



FIG. 33. Schematic figures of the molecular orbits it bond (a) and a bond (b) surrounding 2a 
clusters. 
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FIG. 34. The intrinsic structure of the excited states of 12 Be obtained by VAP calculations 
with the interaction (h). The density distribution of matter, protons and neutrons of the intrinsic 
states are shown at left, middle and right, respectively. 



the the highest spins in the corresponding configurations (0p) 2 (0p)^(0dls) 2 
and (0p)^(0p)^(0dls)l, respectively. 

The 2a + 4n clustering structures are seen in many levels in 12 Be. By an- 
alyzing the single-particle wave functions of the intrinsic structures in similar 
way to 10 Be, we can roughly consider the 0^,1~, and 0^" states as the state 
with 4 neutrons in 7r 2 a 2 , ir^a 1 and 7r 4 orbits surrounding 2a in terms of the 
molecular-orbits ir and a for the valence neutrons. The clustering develops as 
the number of the valence neutrons in the a orbit increases. 

The experimental value for Gamow- Teller transition strength B(GT) from 
the 12 Be ground state to the 12 B ground state is 0.59 |K| deduced from exper- 
imental data of decay. The weaker GT decay from 12 Be than expected with 
the shell-model calculations in the Op-shell model space has been discussed in 
the relations with the possible abnormal structure of the ground state J62|] . 
The experimental B(GT) value is the useful probe to investigate the structure 
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of the 12 Be ground state. One possible reason for the weak decay is considered 
to be that the 12 Be ground state has more intruder admixture of the higher 
shell orbits. However there was no microscopic calculation which can describe 
the weak decay. To evaluate the B(GT) we obtain the wave function for the 
final state by performing VAP calculations of 12 B with (J V ,K') = (1 + ,— 1). 
The calculated B{GT) with the interaction (h) is 0.8 which well agrees to the 
experimental data 0.59. Before the diagonalization for the obtained states, 
the GT strength from the Of state with large 2Huj components is 0.2, and 
that from the OjJ" state with Ohto is 3.0. The reason for the former value is 
that the state with 2Co configurations is forbidden to decay to the n B ground 
state which is almost the normal Ohuj state. The latter value is as much as 
the simple shell-model calculation and is much larger than the experimental 
data. After the diagonalization for the spin-parity eigen states projected from 
the obtained intrinsic states, the 0^ state and the 0^ state are mixed to re- 
distribute the decay strength for the 0hu> state into other states. As a result, 
B(GT) for 12 Be(0^H 12 B(l+) is 0.8 and the one for 12 Be(0^)^ 12 B(l+) is 
2.1. Is it found that the main reason for the weak Gamow- Teller transition 
is that the 12 Be ground state is dominated by the abnormal state with 2hu> 
configurations. On the other hands, the calculations with the interaction (g) 
give much larger B(GT) strength for 12 B(0^)^ 12 B(1 + ) than the experimen- 
tal data because in this case the main component of the ground state of 12 Bc 
is the Qhuj state. 

The abnormal 2p — 2h structure of the ground band of 12 Be may be 
observed in the strength of El transitions. By using the interaction (h) 
which give the abnormal state for the ground band of 12 Be, the E2 transition 
strength in the second band B(E2\ 2^ —> Oj) is as small as the simple AMD 
calculations, however the E2 strength in the ground band B(E2\ 2\ — ► Of) is 
predicted to be almost twice as large as the simple AMD result. 

VI. MECHANISM OF THE DEVELOPMENTS OF CLUSTERING STRUCTURE 

IN BE ISOTOPES 

In the study with the simplest version of AMD and with the VAP calcu- 
lations, it is found that 2a-cluster cores are seen in most of the low energy 
states of Be isotopes heavier than Be. In this section, we try to understand 
the mechanism of the clustering developments in Be isotopes systematically. 
We discuss this problem from two view points. First we take care of the va- 
lence neutrons in the molecular orbits. In the second subsection we give the 
discussion from the other view point of the two-center clusters. 

A. Molecular orbits of valence neutrons 

In the study of Be isotopes with the simplest version of AMD, we have 
already shown the dependence of the development of the clustering structure 
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on the neutron number (Fig. 17). We have discussed the clustering of the 
lowest states of the normal parity states with Ohu> configurations and the 
non- normal parity states with l%uj. In the study of the excited states with 
VAP calculations in the framework of AMD, we have found the rotational 
bands. The interesting point is that the rotaional bands can be classified by 
the number of the neutrons in the sd orbits. 

As already mentioned in the subsection VA, in the case of 10 Be, three 



kinds of the intrinsic states construct the rotational bands. The intrinsic 
states are characterized by the orbits of the two valence neutrons surrounding 
2a cores. The lowest positive parity state in the K = Of band corresponds 
to the state with 2 neutrons in p-orbits. On the other hands, one neutron 
occupies the sd-orbit in the excited states in th K = 1± band. Furthermore 
the states with two neutrons in the sd-orbits appear in the second K = 0^ 
band. The clustering develops larger and larger from K = 0^ to K = V[ and 
to K = 0+. 

According to the VAP calculations, we found these three kinds of the in- 
trinsic states with 2a cluster cores in the other neutron-rich nuclei 11 Be and 
12 Be as well as 10 Be . Namely the most of the excited states in the low- 
energy region are classified into the states with the all the valence neutrons 
in the p-orbits, the states with one neutron in sd-orbit, and the states with 
two neutrons in sd-orbits which surround the 2a cores. Strictly speaking the 
single-particle motion of the valence neutrons surrounding 2a is not neces- 
sarily the ideal independent one but also the correlation between neutrons 
should be important. However, we think it useful to discuss the structure of 
Be isotopes from the viewpoints of molecule-like orbits surrounding 2a be- 
cause the component that neutrons occupy the molecule-orbits is significant 
in the neutron-rich Be isotopes 



Here we improve the figure 17 for the development of the clustering in Be 
isotopes from the view points of the orbits of the valence neutrons surrounding 
2a cores. The Fig.^ presents the schematic figure for the cluster development 
of the lowest states in the given configurations. The lower line is for the state 
with all the valence neutrons surrounding 2a in the p-orbits. The middle and 
upper lines correspond to the state with the one neutron and two neutrons in 



the sci-orbits, respectively. With the help of the figure 35 we can understand 
the development of the clustering in the Be isotopes systematically. 

In the states with the same number of the valence neutrons in sci-orbits, 
the clustering development decreases gradually with the increase of the neu- 
tron number. The mechanism of the clustering development is easily under- 
stood by the idea of molecule-like orbits surrounding 2a cores. In Be isotopes, 
the single particle wave functions of the the valence neutrons surrounding 2a 
cores are analogous to the molecule-like orbits, ir and a orbits, which cor- 
respond to the p-orbit and sd-orbit in terms of the shell-model respectively, 



as shown in the subsection VA. In that sense the additional neutrons in 
each lines occupy the 7r-orbit as the neutron number increases. Therefore it 
seems that the valence neutrons in the 7r-orbits effect as to weaken the spatial 
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FIG. 35. The schematic figure for the development of clustering in Be isotopes. The lines 
correspond to the cluster development of the lowest states in the given configurations (all the 
valence neutrons in the p-orbits, one neutron in the s<i-orbits, and two neutrons in the sd-orbits) 
as the function of the mass number. The intrinsic states of the lowest states well correspond to 
the ones in the rotational bands. 

clustering. 

In the each Be nuclei, as the number of the neutrons in sd orbits increases 
the clustering develops larger and larger. In the case of Be isotopes, the sd 
orbits correspond to the a orbits. It is naturally understood that the valence 
neutrons in a orbits prefer the developed clustering states so as to gain the 
their kinetic energy. That is why the clustering develops as the increase of 
the number of neutrons in the a orbits. 



B. Viewpoint of the di-cluster structure 

As mentioned above, one of the viewpoints for the mechanism of the clus- 
tering development is the explanation with the molecule- like orbits a and n. 
Here we try to understand the clustering in the excited states of 10 Be from 
the other viewpoint of the two-center clustering model. 

The reason for the clustering development in the ordinary nuclei has been 
considered as the system gains the kinetic energy with the development of 
clustering even if the potential energy looses. In order to understand the 
mechanism of clustering in light unstable nuclei, we think it helpful to in- 
vestigate the competition of the kinetic and the potential energies in the 
molecule-like states of 10 Be. 

In the density distribution of VAP results, the structures of the 0^, 1~ 
and O2" states seem to be the two-center clustering structures which consist of 
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6 He+ce (see figure 36). To estimate the dependence of the kinetic and potential 
energies on the degree of the spatial clustering development we represent the 
three kinds of configurations for the 6 He+a system corresponding to the 0^ , 
1~ and O2 states of 10 Be by the simplified AMD wave functions 3>aa/.d(Z) as 
follows. The intrinsic spins of the single-particle wave functions are fixed to 
be up or down for simplicity. For the 6 He+a system with the inter-cluster 
distance d (fm), the centers of single-particle Gaussian wave functions are 
located around two points a\ = (—3d/5y/u,0,0) and 02 = (2d/5y/V,Q,0). 
The three kinds of configurations of the centers for Of , 1~ and Of states are 
shown in the right figures of Fig. ^. We put the centers for p |, p j, n j, 
n 1 at the point a\ and p f" , p J. at a%. The centers for the last 4 neutrons are 
located at the points very close to 02 as 02 ± 5 where enough small 5 is chosen 
so that the angle 9 between S and 02 is 6= tt/2, 7r/4 and corresponding to 
the Of, 1~ and O^ states, respectively. For the parity eigen states projected 
from these three kinds of the simplified AMD wave functions $amd(Z) we 
calculate the expectation values of total, kinetic and potential energies as the 
function of the inter-cluster distance d; 

, m = ((1 ± P)*amd(Z)\H\(1 ± P)$amd(Z)) , . 

{ ' ((1±P)<!> AMD (Z)\(1±P)® AMD (Z)) l ] 

m = ((1 ± P)$amd(Z)\T\(1 ± P)*amd(Z)} (u) 

{ 1 ((1±P)$ AMD (Z)\(1±P)$ AMD (Z)) l } 

/v) = ((1 ± P)*amd(Z)\V\(1 ± P)*amd(Z)) ( ^ 

{ >~ ((1±P)$ AMD (Z)\(1±P)<!> AMD (Z)) ' l ; 

where we omit the spin-orbit and the Coulomb forces for simplicity. In Fig. 



37 we present the total energy, the kinetic energy and the potential energy 
for the three kinds of clustering states 0^, 1~ and 0^~ as the function of the 
distance d between 6 He and a clusters. As shown in the figure for the total 
energy (H), the optimum distances indicate that the clustering structure de- 
velops in the system for the 1~ state and is most remarkable in the state for 
O2 , which is consistent with the present results of VAP calculations. The shift 
of the minimum point of the total energy is understood by the energetically 
advantage of the kinetic part as follows. It is found that the kinetic energies 
in the small d region are sensitive to the configurations, while in the case of 
the potential energies significant differences are not seen in the three config- 
urations. When the cluster approaches each other, the kinetic energy in the 
configuration for the 1~ state becomes larger than 0^ by \7ujJ (about 7 MeV 
for v = 0.17) because all the neutrons in 0^ are in 0s and Op shells. On the 
other hands, in the negative parity state for 1~ one valence neutron must rise 
to the higher sd shell in the small d region. Since in the small d limit the wave 
function is almost same as the harmonic oscillator shell model wave function 
due to the antisymmetrization, the kinetic energy of the sd shell is larger by 
l/2htv than the one of p shell. Thus the system for 1~ loses the kinetic energy 
in the small distance d, as a result, the minimum point of total energy shifts 
to the larger d region than the one in the O^ 1 " state. In the case of the linear 
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FIG. 36. Schematic figures for the intrinsic structure of Of, 1~ and O2" from the view point 
of the two-center cluster model are shown in the left columns of (a), (b) and (c), respectively. 
In the right columns we display three kinds of the configurations of the centers of Gaussians in 
the simplified AMD wavefunctions which indicate the two-center 6 He+a cluster model for the 
excited states 0f(&), l~(b) and 0^~(c) of 10 Be. The black(white) circles correspond to the centers 
of Gaussians of the single-particle wave functions of protons(neutrons). 

configuration for the excited 0^" state, the two valence neutrons occupy the 
sd orbit in the small d limit and the kinetic energy is larger than 0^ by %uj. 
That is why the optimum point d in 0^ is the largest of the three. In other 
words, when the cluster approaches each other, the system feels the repulsive 
force in kinetic part because of the Pauli principle. That is the reason why 
the clustering structure remarkably develops in the 0^~ state. In the analysis 
with this simplified two cluster model, we can conclude that the clustering 
develops so as to gain the kinetic energy. It is compatible with the viewpoint 
of the molecular a orbit. 



VII. SUMMARY 

We studied the structures of the ground and the excited states of light 
unstable nuclei with the theoretical method of AMD. 

Li, Be, B and C isotopes were investigated with the simplest version of 
AMD. We showed the calculated results of binding energies, energy levels, 
radii, magnetic moments, electric moments, and transition strengths com- 
paring with the experimental data. The AMD calculations systematically 
reproduce a lot of experimental data of many nuclei. By analyzing the in- 
trinsic wave functions, we found the drastic changes between shell-model-like 
structures and clustering structures. The theoretical results suggest some new 
features in unstable nuclei such as the clustering structures, the molecular or- 
bits, the opposite deformation between protons and neutrons, the neutron 
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FIG. 37. Total, kinetic and potential energies as functions of the distance between clusters in 
the simplified 6 He+a cluster models for the excited states of 10 Be. The adopted interaction is the 
MV1 force with m = 0.62, and the spin-orbit and the Coulomb forces are omitted. The distance 
d between clusters and the configuration for the Qf , 1~ and 0^ states are defined in the text. The 
solid lines, the dashed lines and the dot-dashed lines correspond to the energies in the system for 
Of , 1~ and 0^ states of 10 Be, respectively. It is $j9md that the minimum point shifts outward from 
Of to 0^ because of the repulsive kinetic force in the small r region. 



skin structures, et al. We discussed these phenomena with the help of the 
experimental informations. 

We adopted VAP calculations to study the excited states of neutron-rich 
Be isotopes. Many excited levels were reproduced and were predicted by 
theory. We showed the theoretical results of the El, E2, and (5 transition 
strength, which agree well with the experimental data. By analyzing the 
structures of excited states in 10 Be and 12 Be, rotational bands K w = 0^, 1 _ , 
0^ are found. Particularly K n = 1~ and 0^ bands have the largely deformed 
states with the molecule-like structures which consist of the developed 2a 
clusters and surrounding neutrons. We extracted the single-particle wave 
functions to discuss the behavior of the valence neutrons, and found that the 
valence neutrons in molecule-like orbits surrounding 2a play important roles 
in the clustering structure. 

In order to understand the mechanism of the clustering development, we 
discussed the clustering structures of Be isotopes from two view points, the 
view point of the molecular orbits of the valence neutrons and the one of the 
two-center cluster model. We concluded that the clustering development in 
neutron-rich nuclei can be understood by the roles of the excess neutrons. 

Although AMD calculations describe many new phenomena seen in un- 
stable nuclei, however, future problems remain to be solved. For example, 
the extremely large radii due to the neutron halo structures are not repro- 
duced because the AMD wave functions written by Gaussians are not enough 
to represent the long tail of the single particle wave functions. It is not ob- 
vious which the effective interaction parameters should be adopted in AMD 
framework. We must more carefully check the properties of the effective in- 
teractions such as the spin-orbit force and the central force. At least when the 
study goes to the heavier nuclear systems, we will soon meet the problems of 
the saturation property. That is to say that it is difficult to reproduce all the 
binding energies of the stable nuclei in the wide mass-number region with a 
set of parameters of the adopted interactions in the present work. Owing to 
the progress of the computational power, it becomes possible to make AMD 
calculations for the heavier unstable nuclei, where we may discover the other 
exotic structures such as the polymers of a clusters and the super deforma- 
tions. 

Finally we stress again that, as far as we know, the present work is the 
first theoretical study which can systematically describe the various clustering 
structures and shell-model-like structures in the ground and the excited states 
of many light unstable nuclei in one framework without assuming any inert 
cores nor the existence of clusters. 

The computational calculations of this work were partially supported by 
the Research Center for Nuclear Physics in Osaka University. We owe the 
calculations to the system in Yukawa Institute for Theoretical Physics in 
Kyoto University, and the super computer systems in High Energy Accelerator 
Research and in the Institute of Physical and Chemical Research. 
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